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4.1 INTRODUCTION 

Concern over the global environmental consequences of fully halogenated chlorofluorocarbons (CFCs) 
has created a need to determine the potential impacts of other halogenated organic compounds on strato- 
spheric ozone and climate. The CFCs, which do not contain an H atom, are not oxidized or photolyzed in 
the troposphere. These compounds are transported into the stratosphere where they decompose and can 
lead to chlorine-catalyzed ozone depletion. The hydrochlorofluorocarbons or hydrofluorocarbons (HCFCs 
or HFCs), in particular those proposed as substitutes for CFCs, contain at least one hydrogen atom in the 
molecule, which confers on these compounds a much greater sensitivity toward oxidation by hydroxyl 
radicals in the troposphere, resulting in much shorter atmospheric lifetimes than CFCs, and consequently 
lower potential for depleting ozone. 

The main objective of this chapter is to review the available information relating to the lifetime of 
those compounds (HCFCs and HFCs) in the troposphere, and to report on up-to-date assessments of the 
potential relative effects of CFCs, HCFCs, HFCs, and halons on stratospheric ozone and global climate 
(through “greenhouse” global warming). 

The lifetimes of the HCFCs and HFCs are determined by their rate of oxidation by hydroxyl (OH) 
radicals in the troposphere. Therefore it was necessary to consider the components of tropospheric 
chemistry affecting the OH radical concentration (i.e., tropospheric ozone, hydrocarbons, and nitrogen 
oxides) as well as the kinetics and degradation mechanism of the halocarbons. Special attention is also 
given to the nature of the products of the degradation of the HCFCs and HFCs being considered as possible 
replacements for CFCs. In parallel with this assessment, an initiative with similar objectives— the Alter- 
native Fluorocarbon Environment Acceptability Study (AFEAS)— was carried out. AFEAS was an in- 
depth examination by over 50 scientists. I he present evaluation has benefited from the availability of this 
material and the conclusions drawn from it. The AFEAS report, Vol. 11, is an Appendix to this report. 


The chapter is divided into three main sections. The first section concerns the oxidation of halocarbons 
in the troposphere. The kinetics and degradation mechanisms of the halocarbons are reviewed and the 
nature and fate of the degradation products are identified. Estimates of current global OH abundances in 
the troposphere are discussed and consequent halocarbon lifetimes are provided. The current picture 
concerning tropospheric ozone and its precursors, as they affect OH, is presented together with a summary 
of recent model predictions of future changes in tropospheric ozone and OH. The second section concerns 
the evaluation of potential effects on stratospheric ozone through the determination of relative (to CFC- 
11) Ozone Depletion Potentials (ODPs), while the third section relates their effects to climate through 
evaluation of relative halocarbon Global Warming Potentials (GWPs) for these compounds. Discussion in 
these sections also emphasizes the uncertainties pertinent to calculations of ODPs and halocarbon GWPs. 


4.2 HALOCARBON OXIDATION IN THE ATMOSPHERE 
4.2.1 Background 


The oxidation of atmospheric trace gases containing carbon, hydrogen, nitrogen, sulphur and halo- 
gens, emitted by a variety of natural and man-made sources, occurs by chemical reactions initiated directly 
or indirectly by solar radiation. These chemical processes yield either non-reactive long-lived products 
such as C0 2 or water vapor, or acidic species such as HNO,, H 2 S0 4 , HF, or HC1 which are removed from 
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the atmosphere by physical processes. The oxidation processes occur mainly in the troposphere, and they 
prevent the build-up of excessive concentrations of these trace substances in the atmosphere. Moreover, 
in the context of the topic of stratospheric ozone depletion, tropospheric oxidation acts as a filter preventing 
the injection into the stratosphere of the source gases for radicals in the NO x , C!O x , and BrO x families, 
which can catalytically destroy ozone. 

The chemical oxidation processes generally, but not exclusively, take place in the gas phase. For 
volatile organic compounds including halogenated hydrocarbons, the major oxidizing species is the hydroxyl 
radical OH, which reacts to abstract hydrogen atoms or adds to unsaturated linkages, thus initiating a free 
radical degradation mechanism. The fully halogenated chlorofluoroalkanes do not contain these reactive 
sites and consequently cannot be degraded in this way. 

The OH radicals are generated photochemically and are maintained at a steady-state concentration 
in the sunlit atmosphere. The average concentration of OH, together with its rate coefficient for reaction 
with a particular compound, determines the atmospheric lifetime of that compound. For a given emission 
rate, the amount of any compound that can build up in the troposphere will be in proportion to its lifetime, 
and it follows that the amount that reaches the stratosphere by transport from the troposphere will be 
heavily influenced by its rate of tropospheric oxidation. Current and future changes in the oxidizing 
efficiency of the troposphere, as measured by the average OH concentration, are clearly of high importance 
in assessing the suitability of hydrochlorofluorocarbons as substitutes for the CFCs. 

The oxidation mechanisms in the troposphere are caused by the presence of ozone, since photolysis 
of ozone in the ultraviolet region leads to the formation of the hydroxyl radical. Moreover, approximately 
10% of the total ozone column is present in the troposphere, and therefore, changes in ozone concentration 
in this altitude range need to be considered as part of the assessment of overall column changes. Tropo- 
spheric ozone also adds significantly to the total greenhouse warming by virtue of its pressure-broadened 
infrared absorption in the atmospheric window, and also by its influence through chemistry on other 
radiative active trace gases (e.g., CH 4 ). 

In recent years there has been a major initiative to understand the budgets and life cycle of tropospheric 
ozone. It is now clear that earlier theories, advocating a purely dynamical description with injection of 
ozone from the stratosphere and removal by surface deposition, are inadequate. In situ chemical production 
and loss of ozone, first suggested in the mid 1970s, clearly plays a role in determining the global tropospheric 
ozone budget and has probably led to increases in 0 3 over large parts of the Northern Hemisphere. The 
reactions controlling ozone production and loss in the troposphere are photochemically initiated and involve 
a variety of trace species including nitrogen oxides, volatile organics, and free radicals in the HO x family 
(OH and HO>). 

There are other oxidizing species present in the troposphere which help to maintain stable composition 
in the atmosphere. These include peroxy radicals, nitrate radicals (N0 3 ), halogen atoms, and ozone. These 
species are generally much less important than OH for tropospheric oxidation. Other oxidizing reactions 
can occur in the aqueous phase in cloud and rain, where hydrogen peroxide (H 2 0 2 ) plays an important role. 
These are not important for oxidation of haiocarbons. 

4.2.2 Tropospheric Chemistry of Haiocarbons 

As discussed in Section 4.2.1 above, most organic compounds including haiocarbons are degraded in 
the troposphere primarily by reaction with the OH radical. The chemical reactivity of haiocarbons is 
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reviewed in this section, considering the dominant atmospheric removal pathways, i.e. gas phase reaction 
with the OH radical and photolysis. Reactions with O, and with the NO, radical are of negligible importance 
as an atmospheric removal pathway. In addition, the specific reaction pathways of the halocarbon degra- 
dation processes prevailing in the troposphere are examined and the principal halogenated degradation 
products are identified. The final fate of these degradation products in the troposphere is discussed by 
considering both their homogeneous gas phase chemistry and their physical removal pathways. The 
reactions of HCFCs and HFCs with Ot'D) atoms are unimportant in the troposphere but may be important 
in producing active chlorine or OH in the stratosphere. 


4.2.2.1 OH Radical Reactions 

The reaction of OH with HFCs and HCFCs involves abstraction of a hydrogen atom to form water 
and a haloalkyl radical (R): 

OH + RH — » H,0 + R 


The rate constants for the reactions of OH with many HFCs and HCFCs have been evaluated as part of 
the AFEAS assessment (see Vol. II, Appendix). Recommendations are given for the five HCFCs and three 
HFCs specified by AFEAS as primary alternatives as well as for all other isomers of C, and C 2 HCFCs 
and HFCs where rate data exist. In addition, recommendations are included for CH,CC1,, CHjCk and 
CH„. The recommended rate constants from this evaluation are given in Table 4.2-1. 

The format used for the presentation of the recommended rate constant data is the same as that used 
by the NASA Panel for Data Evaluation (see DeMore et al., 1987). The rate constant tabulation is given in 
Arrhenius form k(T) = A exp ( - E/RT), and contains the following information: 


1. Reaction stoichiometry. 

2. Arrhenius A factor (in units of cm’ molecule 1 s' 1 ). 

3. Temperature dependence ("activation temperature,” E/R) and associated uncertainty (AE/R). 

4. Rate constant at 298 K (in units of cm’ molecule 1 s' 1 ). 

5. Uncertainty factor at 298 K. 

All of the uncertainties are one standard deviation, lcr. Hence, 95% confidence limits are given by 2cr. The 
uncertainty (lcr) at any temperature can be calculated from the expression: 

f(T) = f(298) exp{A E/R(l/T - 1/298)} 

The rate constants are well defined and significant uncertainties only appear when the Arrhenius expressions 
in Table 4.2.1 are used for extrapolations over a wide temperature range. In particular, at temperatures 
around 277 K, which corresponds to the region of the troposphere where most of the degradation of these 
molecules occurs, the rate coefficients are accurate to within —20%. 


4.2. 2. 2 Photolysis 

Photolysis is a potential atmospheric loss process for those haloalkanes containing multiple Cl and / 
or Br atoms. The absorption spectra are unstructured and continuous, and may extend out to 300 nm, 
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Table 4.2-1. Recommended rate constants and uncertainties for reactions of OH with 
selected HFCs and HCFCs 


Fluorocarbon 


Reaction 

Number 

A a 

E/R ± E/R b 

km* 

f(298) 

OH+CHFCb 

HCFC-21 

1.2(- 12) 

1 100 ± 150 

3.01-14) 

1.1 

OH +CHF.CI 

HCFC-22 

1 .2( — 12) 

1650± 150 

4.71-15) 

1.1 

OH + CHF, 

HFC-23 

7.4(— 13) 

2350 ±400 

2- 1 ( — 16) 

1.5 

OH + CH : CI 2 

30 

5.8( — 12) 

1 100± 150 

1.41-13) 

1.2 

OH + CH ; FCI 

HCFC-31 

3.0(-12) 

1 250 ± 150 

4.51-14) 

1.15 

OH + CH 2 F 2 

HFC-32 

2.51-12) 

1650 ±200 

1.01-14) 

1.2 

OH + CH,F 

HFC-41 

5.41-12) 

1 700 ± 300 

1.8(-14) 

1.2 

OH + CH 4 

50 

2.31-12) 

1700 ±200 

7.71-15) 

1.2 

OH + CHCl,CF, 

HCFC-123 

6.41-13) 

850 ±150 

3.71-14) 

1.2 

OH + CHFCICF, 

HCFC-124 

6.61-13) 

1250 ±300 

1.01-14) 

1.2 

OH + CHF.CF, 

HFC-125 

3.81-13) 

1 500 ± 500 

2.51-15) 

2.0 

OH + CH 2 CICF 2 CI 

HCFC-132b 

3.61-12) 

1600 ±400 

1.71-14) 

2.0 

oh+ch 2 cicf, 

HCFC-133a 

5.21-13) 

1 100 ±300 

1.31-14) 

1.3 

oh+chf 2 chf 2 

HFC-134 

1 

00 

1500 ±500 

5.71-15) 

2.0 

oh+ch 2 fcf, 

HFC-134a 

1.71-12) 

1750 ±300 

4.81-15) 

1.2 

OH + CHjCCI, 

140a 

5.01-12) 

1800 ±300 

1.21-14) 

1.3 

OH + CH,CFCl 2 

HCFC-14lb 

2.71-13) 

1050 ±200 

7.51-15) 

1.3 

OH + CH,CF 2 C1 

HCFC-142b 

9.61-13) 

1650 ± 150 

3.81-15) 

1.2 

oh+ch 2 fchf 2 

HFC-143 

2.81-12) 

1500±500 

1.81-14) 

2.0 

OH + CH,CF, 

HFC-143a 

2.61-13) 

1500 ±500 

1.71-15) 

2.0 

OH + CH 2 FCH 2 F 

HFC-152 

1.71-11) 

1500 ±500 

l.K-13) 

2.0 

oh+ch,chf 2 

HFC-I52a 

1.51-12) 

1 100 ±200 

3.71-14) 

l.l 

oh+ch,ch 2 f 

HFC-161 

1.31-11) 

1 200 ± 300 

2.31-13) 

2.0 


;i Units are cm 3 molecule 's '. 
b Units are K, 


which could result in photodissociation in the troposphere. However, the cross sections are very low in 
the long wavelength region, resulting in very low photolysis rates in the troposphere for fluorochloroalkanes 
and brominated compounds such as CF 3 Br. Only for haloalkanes containing one Br and one or more Cl 
atoms (CF 2 CIBr, CHCFBr, and CH.CIBr) or two or more Br atoms (CF 2 Br 2 , CH 2 Br 2 , CHCIBr 2 , and CHBr 2 ) 
does photolysis in the troposphere become significant. 

Solar photolysis is likely to contribute significantly to the stratospheric destruction of the alternative 
fluorocarbons which have two or more chlorine atoms bonded to the same carbon atom. 

The absorption cross sections for several HFCs and HCFCs have been reviewed as part of the AFEAS 
assessment. Two of the eight HCFCs considered in the review, namely HCFC-123 and HCFC-141b, have 
significant stratospheric photolysis. For these two species there is good agreement among the various 
measurements of the ultraviolet cross sections in the wavelength region which is important for atmospheric 
photodissociation, that is, around 200 nm. There is also good agreement for HCFC-124, HCFC-22, and 
HCFC-I42b; these three species contain only one chlorine atom per molecule and OH reaction is the 
dominant loss process in both the stratosphere and the troposphere. 
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Photodissociation occurs with quantum yields which are expected to be unity, and results in the 
detachment of a Br atom for bromine-containing haloalkanes, or a chlorine atom from HCFCs, e.g., 

CF,ClBr + hv^ CF,CI + Br 
CFiCHCI, + hv CF,CHCI + Cl 

The subsequent reactions of the haloalkyl radicals formed are discussed below. 


4.2.2.3 Reactions of the Haloalkyl (R) Radicals 

Under tropospheric conditions, the only reaction of the haloalkyl radicals, R, is with O,. 

M 

R + O, — » RO, 


For the halomethyl radicals CF 2 , CFC1 2 , and CCI 3 , these additional reactions are in the fall off region 
between second- and third-order kinetics at the temperatures and pressures encountered in the troposphere. 
However, the bimolecular rate constants are estimated to be within a factor of -2 of the limiting high 
pressure rate constants throughout the troposphere, with 

k(R + 0 2 ) s 10 12 cnT molecule ' s~' 
leading to a lifetime of the haloalkyl radicals in the troposphere of <10 6 s. 


4.2.2.4 Reactions of the Haloalkyl Peroxy (R0 2 ) Radicals 

At the expected tropospheric concentrations of NO, N0 2 , the H0 2 radical, and other organic peroxy 
radicals, the dominant reactions of the haloalkyl peroxy radicals will be with NO, NO : , and the H0 2 radical, 
i.e., these reactions proceed by 


RO : + NO-* RO + NO, (a) 

M 

RO, + NO,—* ROONO, (b) 

RO, + HO,-* ROOH + O, (c) 

leading to the formation of haloalkoxy radicals, RO, halogenated peroxy nitrates (ROONO,) and hydro- 
peroxides (ROOH). Based upon kinetic data reported for reactions (a) and (b) of the chlorofluoromethyl 
peroxy radicals (NASA, 1987; IUPAC, 1989), it appears that k a ~1.6 x 10 " (T/300) 1 2 cm 3 molecule " 1 
s 1 and kb — ■ 1 x 10 " cm’ molecule -1 s " 1 for tropospheric conditions. 

No data are available for the reactions of haloalkyl peroxy radicals with the HO, radical. Based upon 
the absolute rate constant data for the reactions of the HO, radical with CH 3 O, and C’HsO,, the rate 
constant for the reactions of alkyl peroxy radicals with HO, is -3.4 x 10 "' 3 e wx ’T cm’ molecule " 1 s', and 
in the absence of experimental data for the haloalkyl peroxy radicals this expression may be assumed to 
apply for reaction (c) above. 
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4.2.2.5 Reactions of the Hydroperoxides (ROOH) 

The gas-phase loss processes of the hydroperoxides are expected to be photolysis and reaction with 
the OH radical. The absorption cross sections for CH 3 OOH have been measured, and photolysis 

ROOH 4 hv — ► RO 4 OH 

is estimated to lead to a lifetime of CH 3 OOH and other ROOH species of a few days in the lower troposphere . 

Reaction of the haloalkyl hydroperoxides, ROOH, with the OH radical will proceed by H atom 
abstraction to reform the RO : radical 


OH 4- ROOH — H 2 0 + RO. 

Abstraction of an H atom from the haloalkyl (R) group, even if possible, is estimated to be of minor or 
negligible importance. Based upon the available kinetic data for (CH 3 ) 3 COOH and CH 3 OOH, the rate 
constant for the reaction of the OH radical with ROOH to form R0 2 4 H0 2 is 1.7 x 10“ 12 e 220 /T cm 3 
molecule 1 s '. This OH radical reaction is then of comparable importance to the photolysis of the ROOH 
species. However, since the estimated lifetime of the ROOH species due to photolysis and reaction with 
the OH radical is ~2-3 days in the lower troposphere, wet deposition and/or incorporation of ROOH into 
cloud, rain and fog water must also be expected to be of importance as a tropospheric sink for these 
hydroperoxides. 


4.2.2.6 Reactions of the Haloalkyl Peroxynitrates (R00N0 2 ) 

The peroxynitrates undergo thermal decomposition and photolysis. Thermal decomposition 

M 

ROONO, R0 2 4 N0 2 

is expected to be the dominant loss process for these compounds in the lower troposphere. Consistent with 
the reverse formation reaction, the decomposition rate constants for the halomethyl peroxynitrates are in 
the fall-off region at the temperatures and pressures of the troposphere, but are within a factor of ~2 of 
the high-pressure rate constants (IUPAC, 1989). For the peroxynitrates CF 2 C100N0 2 , CFC1 2 00N0 2 , and 
CCF00N0 2 , the limiting high-pressure decomposition rate constant is ^ — 1 x I0 l<5 e" ,,00W/T s -1 . 

The absorption cross sections of CFCI 2 00N0 2 and CCl 3 OON0 2 have been measured over the 
wavelength range 210-280 nm (Morel et al., 1980) and the cross sections for these R00N0 2 species are 
generally similar to those for HOON0 2 and CH 3 OON0 2 . Photolysis is expected to proceed by 

ROONO. 4 hv — ► RO, 4 NO, 

ROONO : 4 hv RO 4 NOi 

with an overall quantum yield of unity. The individual quantum yields for these photolysis pathways are 
not known. 
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The tropospheric reactions of the haloalkyl peroxy radicals lead ultimately to the corresponding 
alkoxy radicals, RO. 


4.2. 2. 7 Reactions of the Haloalkoxy (RO) Radicals 

Haloalkoxy radicals undergo three types of reactions under the conditions of the troposphere: 

(1) reaction with 0 : for those radicals containing an a-H atom 

R,R : CHO + Oi —* R,R 2 CO + HO: 

(2) Cl atom elimination for those radicals containing an a-chlorine atom 

R,R : CC10 R,R:CO + Cl 

(3) C-C bond dissociation for s C 2 radicals, 

CX,CX:0 -» CX, + COX: 


where X is an H or halogen atom. 

All of these RO radical reactions yield carbonyl compounds. Experimental studies of haloalkane 
oxidation mechanisms and thermochemical considerations allow the preferred reaction pathway(s) to be 
predicted with a reasonable degree of confidence for each particular RO radical. The tropospheric reactions 
of CF,0 are not well understood, since F atom elimination is too slow to be of any significance. However, 
it is assumed that CF,0 will form COF, under tropospheric conditions. 

The principal carbonyl compounds expected to be formed during the tropospheric degradation reac- 
tions of the haloalkanes are HCOF, HCOCI, HCOBr, COF:, COFCI, COCF, CF,COF, CF,COCI, and 
other CX,COY compounds (X being H and/or halogen and Y a halogen atom). A complete list of the 
products expected from the degradation of the HCFCs and HFCs selected for the AFEAS assessment are 
summarized in Table 4.2-2. 


4.2.2.8 Reactions of Carbonyl Halides, Acetyl Halides, and Halogen-Substituted Acetaldehyde 

Reactions of the carbonyl halides COX: and acetyl halides CX,CfO)Y (X = H and/or halogen; Y = 
halogen) with OH radicals and other reactive tropospheric species are expected to be too slow to be of any 
significance. Furthermore, photolysis is expected to be very slow for these species. 

The halogen-substituted acetaldehydes, CX,CHO (X = F, H, Cl), can, however, react with the OH 
radical 

OH + CX 3 CHO H ; 0 + CX,CO 

to form an acyl radical. As for the acetyl radical (IUPAC, 1989), these CX,CO radicals will rapidly add 0 2 
(with a rate constant of — 10 12 cm’ molecule s - ' under tropospheric conditions) 
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Table 4.2-2. Fluorine-containing products in the atmospheric degradation of selected fluorocarbons 




Atom & 





Compound 

Formula 

Radical 

Carbonyl 

Acid 

Hydroxide 

Nitrate 

HCFC-123 

HCCLCF, 

CF 3 CC1 2 00 

cf,cci 2 o 

CF 3 CCIO 


CF^CCbOOH 

cf,cci,oono 2 



CFiOO 



CF 3 OOH 

cf 3 oono 2 



CF30 



CF 3 OH 

cf 3 ono 2 

HCFC-I4lb 

CCbFCH, 

cci 2 fch 2 oo 

cci 2 fch 2 o 

CCI.FCHO 


CCI 2 FCH 2 OOH 

cci 2 fch 2 oono 2 



cci 2 foo 



cci 2 fooh 

cci 2 foono 2 



cci 2 fo 

CC1FO 






CCI 2 FC(0)00 


CCl 2 FC(0)OOH 

CCI 2 FC(0)OH 


CC1 2 FC(0)00N0 2 

HCFC-142b 

CCIF.CH, 

ccif 2 ch 2 oo 

ccif 2 ch 2 o 

CC1F 2 CH0 


ccif 2 ch 2 ooh 

ccif 2 ch 2 oono 2 



ccif 2 oo 

ccif 2 o 

cf 2 o 


ccif 2 ooh 

ccif 2 oono 2 



CCIF 2 C(0)00 

CCIF 2 C(0)0H 


CClF 2 C(0)OOH 


CCIF 2 C(0)00N0 2 

HCFC-22 

CHC1F 2 

ccif 2 oo 



ccif 2 ooh 

ccif 2 oono 2 



ccif 2 o 

cf 2 o 




HCFC-124 

CHCI- 

CF.CCIFOO 



CF 3 CCIFOOH 

CF3CCIF00N0 2 


FCFi 

CFjCClFO 

CF.00 

CF 3 CFO 


CF 300 H 

CF300N0 2 



CF,0 



CF 30 H 

cf 3 ono 2 

HFC-134a 

CH 2 FCF, 

CF,CHFOO 

CF.CHFO 

CHFO 

CF 3 CFO 


CF 3 CHFOOH 

cf 3 chfoono 2 



CF.OO 



CF 300 H 

cf 3 oono 2 



CFiO 



CF 30 H 

CF30N0 2 



CFO 



CF(0)00H 

CF(0)OON0 2 

HFC-152a 

chf 2 ch, 

ch,cf 2 oo 



CH 3 CF 2 OOH 

ch 3 cf 2 oono 2 



ch,cf 2 o 

chf 2 ch : oo 

cf 2 o 


CHF 2 CH 2 OOH 

chf 2 ch 2 oono 2 



chf 2 ch 2 o 

chf 2 oo 

chf 2 cho 


chf 2 ooh 

chf 2 oono 2 



chf 2 o 

CHFO 






CHF 2 C(0)00 


CHF 2 C(0)OOH 

CHF 2 C(0)0H 

CHF 2 C(0)00N0 2 




CFO 


CF(0)00H 


CF(0)OONO, 

HFC-125 

chf 2 cf, 

cf 3 cf 2 oo 

cf 3 cf 2 o 

cf 2 o 

CFiCFO 


CF 3 CF 2 OOH 

CF3CF 2 OON0 2 



cf,oo 



CF 300 H 

CF300N0 2 



cf,o 



CF 30 H 

cf 3 ono 2 
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M 

CX,CO + O, -»• CX t C(0)0 : 

to form acyl peroxy radicals. These radicals are expected to degrade ultimately to the carbonyl halides and 
carbon dioxide. However, the CX,C(0)0, radicals are expected to form peroxynitrates by reaction with 

NO,: 


CX,C<0)0: + NO, CX,C(0)0,N0, 

The peroxynitrate products are similar to the well-known peroxyacetyl nitrate (PAN) and are expected to 
be thermally stable in the troposphere. In addition, photolysis of these halogenated peroxynitrates is 
expected to be very slow and their residence time due to gas phase degradation could be long. 


4.2.2.9 Physical Loss Processes 


The low gas phase reactivity of many of the halogenated carbonyl compounds produced in the 
degradation of HCFCs and HFCs leads to a potentially important role for physical loss process from the 
atmosphere, i.e., incorporation into rain or clouds, or sea water, with subsequent hydrolysis. Although 
handicapped by the total absence of Henry's law solubility data for any of the compounds of interest and 
the limited availability of relevant kinetic data, an assessment of the rates and mechanisms ol aqueous 
phase removal of the gas phase degradation products has been carried out, as part of the AFEAS assessment. 


The species XXO, HXCO. CH,CXO, CF,OH, CX,OONO„ and ROOH (X = F or Cl, R = halo- 
substituted methyl or acetyl) are all expected to be removed from the atmosphere on time scales limited 
by transport to cloudy regions or the marine boundary layer (i.e., about I month). Some support for t is 
comes from the recent measurements of COCU by Wilson et al. (1988) in the troposphere and lower 
stratosphere which are consistent with a tropospheric physical loss process. Aqueous phase reactions of 
these species result in the formation of chloride, fluoride, and carbon dioxide, as well as formjc acetic^ 
and oxalic acids. The species CX,CXO, CX.CXTXJH, CX,CXT)ONO„ CX,C(0)00NCX, and CX,C(O)0OH 
are also expected to be removed from the atmosphere rapidly, and their aqueous phase reactions resu t in 
the formation of halo-substituted acetates, CX^C(0)0. 


The species CX,C(0)OH are very acidic and, as a result, are highly soluble in cloud water. These 
acids are expected to be rapidly removed from the atmosphere by rainout. However, the aqueous phase 
species CX,C(0)0 are expected to be resistant to chemical degradation. Trichloroacetate can thermally 
decompose on a time scale of 2-10 years to yield carbon dioxide and chloroform. In freshwater, the reaction 
of CCI,C(0)0 , CF,C1C(0)0 may have very long aqueous phase lifetimes. The longest-lived species, 
CF,C(0)0 , could have a lifetime in natural waters as long as several hundred years. Processes which 
could possibly degrade CF n CI, „C(0)0 on shorter time scales than suggested above, but whose rates 
cannot be estimated with any degree of confidence at this time, include oxidation by photochem.cally 
generated valence band holes in semiconductor particles and hydrolysis catalyzed by enzymes in micro- 
organisms and plants; further research aimed at characterizing these processes is needed. 

One possible gas phase degradation product about which very little is known is CF,0N0 2 . This 
compound has never been observed, and may be thermally unstable. If CF,ONO, is thermally stable, then 
it may have a long lifetime toward aqueous phase removal. Henry’s law solubility data and hydrolysis 
kinetics data for CF,0N0 2 are needed before its aqueous phase removal rate can be assessed with any 
degree of confidence. 
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4.2.2.10 Summary of Halocarbon Oxidation Chemistry 

The following main conclusions concerning possible CFC substitutes can be drawn from the discussion 
of halocarbon processes contained in this chapter and in the more detailed AFEAS assessment contained 
in the Appendix to this report. 

Tropospheric reaction with the OH radical is the major and rate determining loss process for the 
MFCs and HCFCs in the atmosphere. The rate coefficients for the OH reactions are well defined at the 
temperatures appropriate for tropospheric reaction. There are virtually no experimental data available 
concerning the subsequent reactions occurring in the atmospheric degradation of these molecules. By 
consideration of data for degradation of alkanes and chloroalkanes, it is possible to postulate the reaction 
mechanism and products formed in the troposphere from HCFCs and HFCs. However, the results are 
subject to large qualitative and quantitative uncertainty, and may even be incorrect. 

A large variety of chlorine- and fluorine-containing intermediate products such as hydroperoxides, 
peroxynitrates, carbonyl halides, aldehydes, and acids can be expected from the degradation of the 
proposed CFC substitutes. Based on the available knowledge of gas phase chemistry, only four of these 
products appear to be potentially significant carriers of chlorine to the stratosphere. These are CCIFO, 
CF,CCIO, CCIF ; CO,NO : , and CC1 2 FC0,N0 3 . However, physical removal processes may reduce this 
potential. In addition, the possibility of pathways and products not predicted by the arguments-by-analogy 
are a cause for concern. 

A large part of the uncertainty of the mechanistic details of the HCFC oxidation arises from an 
insufficient knowledge of the thermal stability and reactivity of halogenated alkoxy radicals. In particular, 
the mechanism of oxidation of the CF,0 radical, which is assumed to produce CF,0, is not known for 
atmospheric conditions and needs further study. Particular attention should also be paid to obtaining data 
on the photochemistry, gas phase reactivity, and solubility of the carbonyl, acetyl, and formyl halides in 
order to assess their removal rates and mechanisms. 

Based on current knowledge, the products identified are unlikely to cause significant changes to the 
effective greenhouse warming potential of the proposed CFC substitutes. This conclusion would be modified 
if long-lived products such as CF,H were formed by unidentified pathways. 

It would be prudent to carry out comprehensive laboratory tests and atmospheric measurements to 
ascertain the validity of the proposed degradation mechanisms for HCFCs and HFCs, before large amounts 
of these substances are released to the environment. 

4.2.3 Hydroxyl Radical Chemistry in the Troposphere 

4.2.3.1 Processes Governing OH Concentration in the Troposphere 

Following Weinstock's (1969) suggestion that reaction with OH radicals provides a major sink for 
atmospheric CO and CH 4 , Levy (1971) was the first to propose that relatively high steady-state concentra- 
tions of OH and other radicals might be present in the troposphere, as a result of photochemical reactions 
involving ozone and other trace gases. Since that time, a comprehensive and relatively consistent theory 
has been developed concerning the sources, chemical cycles, and sinks which control the OH concentration 
in the troposphere. 
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The essential chemistry in the theory of tropospheric OH is illustrated in Figure 4.2-1. The steady- 
state concentration of OH is determined by a balance between production and loss processes for free 
radicals, and the fast interconversion reactions coupling OH with the hydroperoxy radical (HO : ). The main 
production route results from the photolysis ol ozone at \ < 320 nm to produce excited O('D) atoms. The 
O('D) atoms are mainly quenched to the ground state, but a significant fraction reacts with water vapor to 
produce OH: 


O, + hi>(X< 320 nm) -» O('D) + O, 

O('D) + H 2 0 — » 20H 

It follows that the rate of OH production is proportional to the water vapor mixing ratio. In addition, the 
photolysis of H 2 0 2 organic hydroperoxides provide minor secondary sources of OH in the troposphere. 
Because the hydroperoxyl radical, H0 2 , is rapidly converted to OH by reacting with NO and O,, photo- 
chemical reactions leading to formation of H0 2 provide a further, indirect source of OH. The most important 
contribution of this type comes from the UV photolysis of formaldehyde, HCHO, which is present in the 
troposphere as a result of the oxidation of CH 4 and other volatile organic compounds: 

HCHO + hv (A. < 330 nm) — » H + HCO 
HCO + 0 2 -> HO, + CO 

H + 0 2 + M — » HO, + M 

Because CH 4 oxidation is a chain reaction initiated by OH attack (see Figure 4.2-1), the organic oxidation 
chemistry provides an amplification of the primary OH production resulting from ozone photolysis, as 
emphasized by Warneck ( 1975). 



Figure 4.2-1 . Reactions governing concentrations of OH and H0 2 , 
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The loss processes involve combination reactions of the OH and H0 2 radicals with each other and 
also with NO : to form stable, non-radical products, e.g., H 2 0, H 2 0 2 , HN0 3 , and H0 2 N0 2 . The more rapid 
reactions of OH with CO and CH 4 do not lead to overall loss of radicals since H0 2 is formed, which is 
cycled back to OH. Nevertheless, the steady-state concentration of OH is influenced by the concentrations 
of reactive trace gases such as CH 4 , CO, 0 3 , and NO since these determine, through the fast interconversion 
reactions, the partitioning of the total radical pool between OH and H0 2 . 

The possible influence of cloud droplets and aerosol particles in providing heterogeneous loss mech- 
anisms for HO x free radicals has been considered by Warneck (1974). Aerosol loadings in the clean 
troposphere are too low to provide significant sinks; but in the polluted boundary layer and in fogs and 
clouds, heterogeneous loss of H0 2 , which has a longer chemical lifetime than OH, can be an important 
factor in reducing the radical concentration. These factors combined with the dependence of OH on local 
trace gas concentration lead to an expected large local variability in the concentrations of OH. 

Since the dominant primary source of OH is the UV photolysis of ozone, a large seasonal, diurnal, 
and latitudinal variation in the production rate is implied. Moreover, water vapor mixing ratio and cloudiness 
will also influence OH production. These factors were discussed by Warneck (1975) and have subsequently 
been investigated in detail in photochemical models of the troposphere, used to calculate globally averaged 
tropospheric OH concentration (Crutzen and Fishman, 1977; Derwent and Curtis, 1977; Logan et al., 1981). 
The models predict highest OH concentrations in the lower tropical troposphere with very low values in 
polar regions. 

At nighttime, there are no photochemical sources of OH. Nighttime free radical chemistry does occur, 
however, as a result of initiation by attack of nitrate radicals on organic molecules, but OH steady-state 
concentrations are much lower than in the sunlit atmosphere. 


4.2.3.2 Validation of OH Photochemistry 

In contrast to the well-developed theory of tropospheric OH, experimental measurements of OH 
concentrations and verification of the theory have provided a challenge which has not yet been satisfactorily 
overcome. 

Three techniques for the time-resolved measurement of local OH concentrations have been developed 
over a number of years, these being the l4 C-tracer method, laser-induced fluorescence (LIF), and long-path 
absorption spectroscopy (LPAS). Both the l4 C-tracer method and the LIF technique are subject to a number 
of potentially severe interferences, and the results of even the most recent measurements using these 
techniques must be treated with caution. To date the UV LPAS technique appears to give the most reliable 
measurements. The most recently developed high resolution instrument (Dorn et al., 1988; Platt et al., 
1988) overcomes interference problems by detecting a spectral range covering several OH absorption 
features. This method gives absolute OH concentrations and does not require calibration. The only 
information needed for data analysis is the oscillator strength for OH which is known to be better than 5%. 

The daytime OH measurements published over the last four years using all three techniques vary 
between approximately 0.4 and 9 x I0 h OH radicals per cm\ which is broadly consistent with current 
model predictions. Validation of a time-dependent model, however, requires not only the measurement of 
OH itself, but also the simultaneous determination of all input parameters that control the local OH 
chemistry. In the recent LPAS measurements, sufficient supporting data were obtained to allow full model 
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calculation of the expected local OH concentrations (Perner et al., 1987). Comparison of these recent 
experimental results with the model calculations leads to the following conclusions: 

1. in the presence of high NO; concentrations (>2 ppbv), the steady-state OH concentration is 
controlled only by the primary source term (ozone photolysis and subsequent reaction of O('D) 
atoms with water vapor) and one main loss reaction (OH + N0 2 forming nitric acid). Thus the 
calculation of OH concentrations needs only measurements of O,, J(O'D), H : 0, and NOj, with 
the most critical being the correct determination of J(O'D). In this case the measurements agree 
with the calculated OH concentrations to within —50%. 

2. In relatively clean air (N0 2 < 1 ppbv), the OH concentration calculated using a detailed model 
with the supporting measurements as input parameters exceeded the measured concentrations 
by about a factor of two. This implies an inadequate description of the HO x losses in the model. 

Further development of techniques for OH measurements arc clearly necessary and confidence in 
the measurements would be aided by a successful intercomparison of more than one system. Also, there 
is a need for measurements at different characteristic locations, in particular at high altitudes remote from 
local terrestrial sources of gases and particles. 

4.2.3.3 Indirect Determination of Global Hydroxyl Radical 

Indirect methods for determining the OH abundance from trace gas budget considerations have 
provided the best test for global models of OH distribution. In principle, any molecule can be used provided 
that the following conditions are satisfied. 

• Its removal occurs by reaction with OH alone and the rate of this removal is accurately known 
both as a function of temperature and pressure. 

• Its atmospheric abundance is accurately known. 

• The source strength and distribution in time and space is accurately known. 

While many chemicals (e.g., CO, non-methane hydrocarbons [NMHC]) are nearly exclusively removed 
by OH radicals, it has not been possible to use them because their sources are diverse and often unknown. 
However, in two cases the above conditions are met to varying degrees: man-made halocarbons and 
naturally occurring l4 CO . In recent decades synthetic halogenated chemicals have been injected into the 
environment in such significant amounts that a measurable global background is present. Because of the 
exclusive man-made source of halocarbons, the major uncertainties associated with the source term can, 
in principle, be eliminated. A select group of these halocarbons are also removed from the atmosphere 
nearly exclusively by reaction with OH. Similarly, l4 CO has a well-defined cosmic ray source and its 
atmospheric abundance can be related to OH. In the following paragraphs previous attempts to evaluate 
OH abundance in this way are summarized and a recent assessment of the global OH concentrations and 
resultant halocarbon lifetimes conducted within AFEAS is presented. 

Methyl chloroform and the hydroxyl radical 

A first demonstration of this technique was presented by Singh (1977a, b) and Lovelock (1977), who 
selected methyl chloroform (CH,CCld as the currently most suitable OH-tracer fulfilling the above criteria. 
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Analysis of CHiCCP data have been performed by numerous investigators during the last decade. In 
most cases, box models of various kinds have been used to estimate its lifetime. Table 4.2-3 summarizes 
the global average lifetime of CHiCCI* derived in a variety of published studies. The first emission inventory 
reported by Neely and Plonka (1978) has been widely used. Subsequent emissions data have been estimated 
from production figures assuming that the production-to-emission ratio of Neely and Plonka ( 1978) has not 
changed significantly. Only 3 to 4% of the emissions are thought to occur in the Southern Hemisphere. 
Prinn et al. (1987) estimate that the CHtCCIt global emission data are probably known to be —5%. A recent 
analysis of methyl chloroform emissions by Midgely (1989) confirms this. 

A bulk of the uncertainty in CH^CCU lifetime reported in Table 4.2-3 comes from the determination 
of the atmospheric burden of CH<CCP. Over the years the estimates of the global average CH,CCI., lifetime 
cover the range 6-12 years. The most recent ALE/GAGE data after correction (see Table 4.2-3) result in 
a mean CH^CCIi lifetime of 6 to 7 years. 

The CHjCClj mean lifetime estimate can be compared to that predicted by a global model with a well- 
defined OH field. It is convenient to define an equivalent average OH concentration that is consistent with 
the CHiCCl, global lifetime. Such an estimate of OH must take into account the fact that OH + CHjCCl* 
reaction rate is a function of temperature and also must assume that other removal processes are negligible. 
With these assumptions a mean OH concentration of 5-7 x 10 s molecule cm ; fits the data best. The 
hemispheric asymmetry in the CH,CCI, lifetime (Table 4.2-3) also implies that the OH levels in the SH are 
somewhat higher than the NH values. 


Table 4.2-3. Atmospheric lifetime of methyl chloroform and estimates of equivalent OH 
concentrations 


Global average 
CH 3 CCI 3 
Lifetime 
(years) 

Mean OH 

(10 s molecule cm 3 )“ 

Comments 

Reference 

7 ± 1 

3-6 


Singh (1977a) 

5-10 

5- 10 


Lovelock (1977) 

8-11 

4 

Tnh 1 5 Tsh 

Singh (1977b) 

3 

10 

Tnh ~ 3 T S h 

Neely and Plonka (1978) 

8 

6 


McConnell and Schiff (1978) 

7 ± 1 

— 


Makide and Rowland (1981) 

5 ±2 

7±2 


Derwent and Eggleton (1981) 

9± 2 

6 ± 1 


Singh et al. (1983) 

iorv ,) b 

5 ±2 


Prinn et al. (1983) 

6 ± 1 .5 C 

8 


Khalil and Rasmussen (1984) 

6± P 

7 ±2 


Fraser et al. (1986) 

12 h 

— 


Golombek and Prinn (1986) 

6 ± P 

8 ± 1 

T"nh ~ P 15 Tsh 

Prinn et al. (1987) 

6.5 ± 1 

— 

Tnh ~~ 1 -2 T sm 

Blake and Rowland (1988) 


a Using model results it is reasonable to assume that the OH abundance between 0-30° latitude is nearly twice as large as 
between 30-90° latitude. Thus, a global OH average (molecule cm A ) of 7.5 x 10‘’ would correspond to OH values of 10 x 10 <J 
and 5 x 10 s at 0-30° and 30-90° latitudes, respectively. 
h Uncorrected ALE/GAGE CH 3 CCI ; , data. 

^Corrected ALE/GAGE CH 3 CCI 3 data. All atmospheric measurements multiplied by 0 8 The actual derived lifetime is 6.3 years. 
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Other halocarbons have been used in a similar manner to estimate OH. for example, dichloromethane, 
1,2 dichloroethane and tetrachloroethene (Singh et ah, 1983). A two-box model shows that removal rates 
of these molecules are consistent with average OH concentrations ot 4 to 5 x 10 s molecules cm '. Because 
of their relatively short lifetimes (several months) compared to CH,CC1,. these chemicals could potentially 
reveal OH latitudinal and seasonal gradients with greater sensitivity. However, their source strengths are 
not as well defined as that of CH,CCI,. In addition, these analyses are sensitive to the seasonal variations 
in emissions which are poorly known. 


14 CO and the hydroxyl radical 

The first attempts to calculate the global mean tropospheric hydroxyl concentration from isotopic 
distribution in atmospheric carbon monoxide were made by Weinstock (1969) and by Weinstock and Niki 
(1972). They used the three available l4 CO measurements by McKay et al. (1963), an estimate of the global 
source strength for M CO from the well-known cosmic ray bombardment of atmospheric nitrogen molecules, 
and the OH + CO rate coefficient to derive an estimate for the atmospheric turnover time for i: CO of the 
order of 1 month. Furthermore, they suggested that OH radicals were responsible for the CO removal and 
obtained an estimate of the global mean OH abundance. 

The ,4 CO concentrations in the troposphere are extremely small with a winter maximum of 20 molecule 
cm 4 and a summer minimum of about 10 molecule cm '. Volz et al. (1981) improved the methodology 
based on four refinements, viz: 

• additional measurements of l4 CO in the lower troposphere, 

• evaluated chemical kinetic data for the OH + CO reaction, and 

• a global two-dimensional time-dependent model to investigate the coupled CH 4 -H 2 -CO-NO x -0, 
life cycles, replacing the box model approach. 

These analyses led Volz et al. (1981) to obtain a mean tropospheric OH concentration of 6.5 4, - 2 x 10' 
molecule cnrl No additional data on l4 CO have been collected in the intervening years and no Southern 
Hemispheric data are available. Derwent and Volz-Thomas (1989) in a re-evaluation of this technique 
conclude that the average OH derived by this method remains unchanged. 


4.2.3.4 Recent Assessment of Tropospheric OH Abundance and Halocarbon Lifetimes 

The above summary shows that an overall volume-averaged OH concentration of 6 (±2) x 10 s 
molecule cm ' best represents all of the available information on the CH,CCI, and l4 CO budgets. However, 
the usefulness of a global mean OH concentration for determination of overall reaction rates and lifetimes 
of HCFCs derived in this fashion is questionable. This arises because of the non-uniformity of the OH 
abundance and of the distribution of the molecules with which it reacts, and also because of the different 
temperature dependencies of the reaction rates in the troposphere. In the recent AFEAS assessment these 
factors were considered at some depth. Prather (1989) used 3-dimensional tropospheric OH fields that were 
calculated from a climatology of sunlight, temperature, and trace gas mixing ratios. These OH fields were 
then used to study the methyl chloroform budget by calculating the integrated loss of this molecule in a 3- 
D chemical transport model of the troposphere and to test the accuracy of scaling the HCFC lifetimes to 
an assumed methyl chloroform lifetime. 
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In the AFE1AS assessment, the lifetimes of HCFCs and HFCs were determined by three separate 
approaches: 

1. 2-D chemical transport model with semi-empirical fit to l4 CO: 

2. Photochemical calculation of 3-D OFI fields and integrated loss; 

3. Scaling of the inferred CHiCCh lifetime by rate coefficients. 

Resulting lifetimes from all three independent approaches generally agree within 15%, as shown in 
Table 4.2-4. The integrated losses calculated from the global OH fields in the models ( 1 and 2) are constrained 
by modeling of the observations and budgets for ,4 CO and CH ,CCI , respectively. Method (3) may be 
expressed simply as 

lifetime(HCFC) = 6.3 yr x k(CH,CCl, at T = 277 K) / k(HCFC at T 277 K), 

where the current estimate of the lifetime for methyl chloroform (6.3 yr) is based on the ALE/GAGE 
analysis (Prinn et al., 1987). Some of the errors associated with this scaling and the choice of the temperature 
of 277 K for scaling the reaction rates have been tested with the 3-D OH fields from method (2). Based on 
the sensitivity studies from methods ( 1) and (2), method (3) should be reliable for calculating HCFC lifetimes 
in the range I to 30 years. Method (2) shows that the middle tropical troposphere (2-6 km) dominates the 
atmospheric loss and would be an important region in which to make observations of OH. 


Table 4.2-4. Atmospheric lifetimes for HCFCs and HFCs, based on Derwent and Volz-Thomas (1989) 
and Prather (1989) 


Lifetime (yr) for method 8 


HCFC 

k(cm 3 molec ' s') 

(D 

(2) 

( 3 ) 

CH,CCI, 

5.0 x 10 

12 exp( — 1800/T) 

5.0 

5.4 

6.3 

(range) 



(3-7) 

(4-7) 

(5.4 — 7.5) 

CH,F 

5.4 x 10 

12 exp( - 1700/T) 

3.3 

3.8 

4.1 

ch 2 f 2 

2.5 x 10 

12 exp( - 1650/T) 

6.0 

6.8 

7.3 

CHF, 

7.4 x 10 

13 exp( - 2530/T) 

635.0 

289.0 

310.0 

CH : FCI 

3.0 x 10 

12 cxp( “ 1250/T) 

1.26 

1.33 

1.44 

CHFCF 

l .2 x 10 

12 exp( - 1 100/T) 

1.80 

1.89 

2.10 

CHF 2 C (22) 

1 .2 x 10 

12 exp( - 1650/T) 

13.0 

14.2 

15.3 

CH,CH 2 F 

1.3 x 10 

" exp( - 1200/T) 

0.31 

0.25 

0.28 

CH.FCH : F 

1.7 x 10 

11 exp( - 1500/T) 

0.60 

0.58 

0.63 

CH,CHF : (l52a) 

1 .5 x 10 

12 exp( - 1 100/T) 

1.46 

1.53 

1.68 

CH.FCHF, 

2.8 x 10 

12 exp( - 1 500/T) 

3.2 

3.5 

3.8 

CH.,CF, 

2.6 x 10 

" exp( - 1500/T) 

40.0 

38.0 

41.0 

CHF : CHF : 

8.7 x 10 

'■’exp(- 1500/T) 

10.4 

11.4 

12.3 

CH.FCF, (134a) 

1.7 x 10 

12 exp( - 1750/T) 

13.1 

14.4 

15.5 

CHF.CF, (125) 

3.8 x 10 

11 expo 1500/T) 

25.9 

26.1 

28.1 

CH,CFCI : (141b) 

2.7 x 10 

" exp( - 1050/T) 

6.7 

6.7 

7.8 

CH,CF,ci(142b) 

9.6 x 10 

exp(- 1650/T) 

16.6 

17.8 

19.1 

CH,CICF : CI 

3.6 x 10 

12 expo 1600/T) 

3.5 

4.0 

4.2 

CH.CICF, 

5.2 x 10 

" exp(-l 100/T) 

4.1 

4.4 

4.8 

CHCFCF, (123) 

6.4 x 10 

” exp( - 850/T) 

1.40 

1.42 

1.59 

CHFCICF, (124) 

6.6 x 10 

13 exp(- 1250/T) 

5.5 

6.0 

6.6 


‘Lifetimes from Method (1) do not include stratospheric loss; those from Method (2) include small additional stratospheric loss. 
Method (3) is based on scaling the methyl chloroform lifetime of 6.3 yr (Prinn et al. 1987) by the ratio of the rate coefficients at 
2 77 K 
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Estimated uncertainties in the HCFC and HFC lifetimes between I and 30 years are ± 50% for (1) 
and ± 40% for (3) and (2). These ranges include the uncertainty in the rate constants for the OH reactions. 
Global OH values that give lifetimes outside of these ranges of uncertainty are inconsistent with detailed 
analyses of the observed distributions for l4 CO and CH,CCI,. The expected spatial and seasonal variations 
in the global distribution of HCFCs with lifetimes of I to 30 years have been examined with methods (1) 
and (2) and found to have insignificant effect on the calculated lifetimes. Larger uncertainties apply to gases 
with lifetimes shorter than one year; however, for these species our concern is for destruction on a regional 
scale rather than global accumulation. 

Clearly any changes in the concentrations of ozone and other trace gases in the troposphere, which 
affect OH abundance, will have an impact on HCFC lifetimes and in the following section we review the 
relevant tropospheric chemistry and current perception of future atmospheric behavior in this context. 

4.2.4 Tropospheric Chemistry Influencing OH and Ozone 

4.2.4.1 Processes Controlling Tropospheric Ozone 

Tropospheric ozone may be produced by in situ chemistry (Crutzen, 1973; Chameides and Walker, 
1973; Fishman and Crutzen, 1978) or by transfer from the stratosphere, where O, is generated by the 
photodissociation of molecular oxygen at altitudes above 25 km, followed by combination of the ground 
state oxygen atoms with CL: 


O: + hv (X< 245 nm) -* 20('P) 

O('P) + O, + M — » O., + M 

The in situ source of tropospheric ozone is the reaction of HO, or organic peroxy radicals with NO 
followed by the photodissociation of the nitrogen dioxide produced and the O('P) + O, reaction. 

HO: + NO —* NO, + OH 
NO, + hv (X < 400 nm) -> O(T) + NO 

Although NO and NO, concentrations are very low throughout most of the troposphere, the effect of 
this reaction sequence can provide a significant O, source. 

Ozone is also consumed in the troposphere as a result of pholodissociation in the ultraviolet region. 
Most of the excited O('D) atoms are quenched by N, and 0 2 to the ground state, OCP), which reform 
ozone, but that fraction which reacts with other trace molecules (e.g., with water to produce OH), represents 
a net loss of ozone. Photolysis of ozone in the near U V (Huggins bands) and visible (Chappuis bands) does 
not lead to net ozone loss. There is also an in situ loss of ozone through reaction with NO : (to give NOT, 
HO: and unsaturated hydrocarbons. Model calculations (e.g.. Levy et al., 1985) point to an approximate 
balance between in situ sources and sinks for tropospheric ozone, averaged over the globe. This is consistent 
with estimates of the magnitude of the surface sink for ozone, which averaged over the globe, is approxi- 
mately equal to the stratospheric injection. However, analysis of observations and regional budgets (Logan, 
1985; Bojkov, 1986; Penkett, 1988) show a clear indication that Northern Hemisphere tropospheric ozone 
is increasing. Changes in the seasonal modulation of ozone concentration, and in the man-made emissions 
and concentration of trace gases which affect ozone, suggest that the increase is due to a shift in the balance 
between in situ production and loss (Penkett, 1988). This balance essentially depends on the amount of 
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NO, present in the background atmosphere together with the concentration of peroxy radicals derived 
from oxidation of CO, CH 4 and non-methane hydrocarbons, which controls the flux through the reactions 
in which NO is oxidized to N0 2 . 

RO : + NO — * NO : + RO (R = H or organic radical) 

The amount of tropospheric ozone is clearly central to the problem of the oxidizing efficiency of the 
troposphere, since 0 3 photolysis is the primary source of OH radicals as well as being an oxidizing species 
itself. It follows that the future trends of oxidizing capacity will be tied to the future tropospheric burden 
of CO, CH 4 , non-methane hydrocarbons and other organics, as well as nitrogen oxides. We now summarize 
the current picture of the distribution of ozone precursors in the troposphere. The distribution of ozone 
itself in the troposphere is covered in the assessment of global ozone concentrations in Chapter 2. 


4. 2. 4. 2 Tropospheric Distribution of Reactive Nitrogen Species 

The tropospheric lifetime of NO, (NO + NO : ) is rather short and consequently it is not possible to 
define an average, representative NO, mixing ratio on a hemispheric or global scale. 

The NO, levels in urban and industrial areas are determined largely by anthropogenic sources. This 
leads to near surface NO, mixing ratios in rural areas of the eastern United States and Western Europe 
typically in the range between 1 ppbv and 10 ppbv. In less populated and coastal regions, the NO, levels 
depend on the prevailing meteorology and the proximity and distribution of urban sources. In these locations 
the NO, mixing ratios typically range between 0.1 and 1 ppbv. 

Surface measurements at remote locations reflect the influence of natural sources, such as soil 
emissions and lightning, which are principally terrestrial sources that are seasonal. These sources may 
provide surface concentrations of NO, that typically range between 0.02 ppbv and 0.1 ppbv (Fehsenfeld et 
al., 1988). It should be noted that at most locations on land, particularly in the mid-latitudes, the NO, levels 
in ambient air are influenced by anthropogenic activities. In remote maritime air and in the polar regions 
that are not influenced by anthropogenic activities, the NO, concentrations are exceedingly small, typically 
0.001 ppbv to 0.01 ppbv, and are associated with the downward mixing of NO, from NO, reservoirs in the 
upper troposphere. 

The free tropospheric burden of NO, is also strongly influenced by anthropogenic NO,, particularly 
from combustion sources in the Northern Hemisphere. Variability in the atmospheric transport and pho- 
tochemical lifetimes, together with the contribution of aircraft emissions of NO, and natural sources such 
as lightning and stratospheric subsidence, increases the variability of NO, concentration in the free tro- 
posphere. Here the NO, mixing ratios may vary from 0.02 ppbv in remote regions to 5 ppbv over populated 
areas. In general, throughout the free troposphere the NO, mixing ratio increases with height, as illustrated 
by the NO distribution shown in Figure 4.2-2. 

The distribution of NO y (NO y : sum of all oxidized nitrogen species except N 2 0 = NO + N0 2 + 
N0 3 + 2 x N 2 0 5 + HNO,-f PAN + H0 2 N0 2 + organic nitrates) is very similar to that of the source 
compound NO,. Near the sources the majority of NO y is in the form of NO,. However, the NO, will be 
rapidly converted to PAN, HN0 3 , or other NO y compounds. For example, during periods of maximum 
insolation the photochemical lifetime of NO, is less than one day. Thus, in the remote free troposphere 
NO, accounts for only a small fraction (approximately 10%) of the NO y . For this reason, NO y is a more 
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MIXING RATIO (ppt) 

Figure 4.2-2. Schematic presentation of the measured vertical distribution of NO and NO y over midlatitude 
oceanic and coastal regions. 


conserved compound and the variation in the concentration of NO y will be less than that of its short-lived 
precursor, NO*. 

In the remote troposphere the main fraction of the total NO y often consists of HNO, and PAN (Fahey 
et al., 1986). HNO, can be quite effectively removed from the atmosphere by wet and dry deposition, 
whereas PAN is reasonably quickly thermally decomposed back into peroxyacetyl radicals and NO ; at 
temperatures representative of the atmospheric boundary layer at low latitudes and at mid-latitudes in 
summer. Due to its stability at lower temperature, sometimes PAN levels of several hundred pptv can be 
observed in the free troposphere at mid- and high northern latitudes (Singh et al., 1986). Thus PAN and 
possibly other organic nitrates can serve as a potential reservoir and transport medium for the rather short- 
lived NO x . However, up to now no quantitative estimate of the possible magnitude of such an effect has 
been made. 

There is evidence from analysis of ice cores from Greenland and the Alps of a substantial increase in 
the concentration of nitrate ions, formed from NO ; oxidation, over a similar period (1895-1978) (Neftel et 
al., 1985). This points to an increasing concentration of NO x in parts of the Northern Hemisphere tropo- 
sphere. Because of their short lifetimes, any increases in nitrogen oxides are expected to be confined to 
the continental source areas and nearby regions. The influence of nitrogen oxides may be spreading, 
however, due to a widening of source areas associated with increasing economic activity in the developing 
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world and due to reservoir species such as peroxyacetyl nitrate (PAN), which are formed in NMHC-NO x 
interactions and which have relatively long atmospheric lifetimes and can dissociate to produce N0 2 far 
from source regions. 

4.2.4.3 Tropospheric Distribution of Hydrocarbons and CO 

The most important and abundant atmospheric hydrocarbon is methane. Its global distribution is well 
established and it exhibits only a relatively small variability in the background troposphere. The Northern 
Hemispheric CH 4 mixing ratios are about 1.7 ppm, the Southern Hemispheric levels are slightly lower, 
about 1 .6 ppm. 

CO has a considerably lower tropospheric abundance, but due to its higher reactivity towards OH 
radicals it is very important for the oxidizing efficiency of the atmosphere. The tropospheric distribution 
of CO is reasonably well known and there is a significant interhemispheric gradient. CO mixing ratios in 
the Southern Hemisphere are around 50-70 ppbv, and roughly a factor 2-3 higher in the Northern Hemi- 
sphere. A distinct seasonal variation in CO levels in both hemispheres has been identified. However, global 
coverage of CO distribution is still inadequate. 

The current picture of the distribution of non-methane hydrocarbons (NMHC) in the troposphere is 
still far from complete due in part to the limited number of investigations, and there is large variability in 
the background concentrations. This variability results from the relatively short average atmospheric 
residence times, which range from a few' hours to a few months. Also, the different NMHCs have various 
sources with different geographical distributions. The most important sources are listed in Table 4.2-5. 

Large-scale distributions of NMHC in the remote troposphere have been published by several inves- 
tigators (Singh and Kasting, 1988; Rudolph, 1988). Most of these studies focused on the latitudinal variability 
of NMHC, but there are also a few studies which present data on the vertical distribution of NMHC. 


Table 4.2-5. Estimates of global hydrocarbon emissions into the atmosphere 3 


Source 

c 2 h 6 

c 3 h 

C 4 H,o 

c 2 h 2 

C 6 H 6 

c 2 h 4 

c 3 h 6 

c s h 12 

c 7 h 8 

CaH 10 

CsH 8 

CioH 16 

Engine exhaust 

+ 


+ 


+ 


+ + 

+ + 

+ + 

+ 

+ 

+ + 

+ + 



Evaporation losses 





+ 

+ 


+ 



+ + 

+ 

+ 



Natural gas leakage 

+ 

-1- 

+ 

+ 

+ 

+ 










Oil + coal burning 
















Chemical industry 

+ 








+ 



+ 




Solvent use 












+ + 

+ 



Biomass burning 

+ 

+ 

+ 

+ 



+ 

+ 

+ + 

-1- + 


+ 

+ 

+ + + 

+ + + 

Foliage emissions 

+ 








+ 





+ 4- + 

-b + + 

Microbial production 

+ 








+ + 

+ 






Ocean emissions 

+ 

+ 

+ 

+ 

+ 

+ 



+ + + 

+ + + 

+ + 





a Data from Ehhalt and Rudolph, 1984. 
Note: + = Moderate source 

+ + = Strong source 

+ + + = Very strong source 
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The distribution of the longer-lived NMHC (t > I week) show some systematic features which can 
be ascribed to a representative latitudinal profile. On the average the highest mixing ratios are observed at 
mid- to high northern latitudes, roughly 1-3 ppbv of ethane, 0.2-0. 6 ppbv of propane and n-butane, 0.1- 
0.5 ppbv of acetylene, 0.05-0.25 ppbv of benzene, and 0.02-0.1 ppbv of i-butane. In general, all these 
species show a considerable decrease towards lower latitudes and this gradient is more pronounced for the 
shorter-lived of these NMHC. This reflects both the source distribution with high emissions in the indus- 
trialized zone of the Northern Hemisphere and the faster removal at low latitudes due to the higher OH 
radical concentrations at tropical latitudes. Average Southern Hemisphere mixing ratios are roughly 0.15- 
0.4 ppbv ethane, 0.03-0.1 ppbv propane, 0.01-0.05 ppbv n-butane, 0.01-0.05 ppbv acetylene, 0.01-0.03 
ppbv benzene, and 0.005-0.03 ppbv i-butane. There is a clear seasonal cycle in the troposphere mixing 
ratios of the light aliphatic hydrocarbons, with substantially higher values in the winter months. The cycle 
is consistent with predominant removal of these hydrocarbons by reaction with OH. 

There seems to be a slight but systematic decrease from low southern latitudes towards mid southern 
latitudes. For ethane and probably the other hydrocarbons, there seems to be a significant and systematic 
seasonal cycle. In both hemispheres the minima are in the respective late summer, maxima in late winter. 

For more reactive NMHC with average atmospheric residence times of less than one week, it is no 
longer justified to consider their distributions as being systematic with latitude. These atmospheric con- 
centrations are determined by emissions, removal reactions, and transport on a local or regional scale (up 
to several hundred km). The most important NMHC in this category (I week > t > 0.5 days) are C;— C 5 
alkenes and, to a lesser extent, < C 5 alkanes and alkylbenzenes. Since the sources of C:-Q alkenes include 
biomass burning, emissions from vegetation, and emissions from oceans, substantial concentrations in the 
range from 0. 1 to a few ppbv of ethene and propane have been observed over nonindustrialized continental 
regions and the oceans. 

The extremely reactive terpenes and isoprene in general show a very strong decrease with increasing 
altitude and can generally only be observed near their sources within the atmospheric boundary layer. 
Over areas with dense vegetation, isoprene mixing ratios around 0.5-2 ppbv and monoterpene mixing 
ratios approaching I ppbv are frequently observed (Rasmussen and Khalil, 1988; Zimmerman et al., 1988). 
The emissions of many biogenic NMHC strongly depend on temperature, relative humidity, and light 
intensity. The types of emitted compounds also depend on the type of vegetation, e.g., isoprene is primarily 
emitted from deciduous plants whereas coniferous trees mainly act as terpene sources. The estimated 
global emission rate of isoprene is about 480 x 10' 2 g yr 1 (Rasmussen and Khalil, 1988). Comparable 
emission rates have been estimated for the other terpenoid compounds. 

The global source strength of all NMHCs, but not including isoprene and the terpenes, is estimated 
to be about 100-130 x 10 12 g yr '. Since these NMHC emission estimates are highly uncertain and probably 
underestimated, we cannot yet make any quantitative estimates on the contribution of NMHC to the 
tropospheric O, and OH concentrations. Also, the atmospheric oxidation of NMHC leads to the in situ 
production of CO, which has a longer residence time than the NMHC and could extend their influence 
over a wider spatial regime. A comparison with the global CH 4 emission rate of about 400 x 10 12 g yr~' 
shows that the possible contribution of NMHC (not including isoprene and terpenes) can be quite substan- 
tial. The influence of isoprene and the terpenes on O, and OH may also be very substantial but, as yet, this 
effect has not been quantified. 

The available measurements of non-methane hydrocarbons in the troposphere do not allow a system- 
atic long-term trend to be recognized. Moreover, the observed variability of NMHC will prevent easy 
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recognition of any global trend in the near future. It appears more promising to estimate future (and past) 
changes of atmospheric NMHC concentrations from possible changes in their main sources. The fossil fuel 
activities have led to large increases in the emissions of non-methane hydrocarbons (NMHCs), but as yet 
these are on a smaller scale than methane itself and again the magnitude of the sources is in dispute. The 
question of the influence of the breakdown of halogenated organic molecules such as the HCFCs and HFCs 
on tropospheric ozone production has been reviewed in the AFEAS assessment (Niki, 1989). It is concluded 
that these molecules are such a small fraction of the total budget of organic compounds oxidized that their 
effect will be trivial. 

4.2.5 Model Evaluation of Trends in Tropospheric Ozone and OH 

As shown in Chapter 2, there is now observational evidence that changes are taking place in the 
concentration of tropospheric ozone and its precursor gases such as CH 4 , CO, and possibly NO*. The 
hydroxyl radical concentration is controlled by fast gas phase reactions. Simple model analysis of the 
photochemistry (Cox and Derwent, 1981 ; Liu, 1988) shows that increases in CH 4 and CO tend to decrease 
OH concentration whereas increases in NO x and O^ tend to increase OH. Changes in the solar UV due to 
depletion of stratospheric ozone could potentially lead to an increase in the oxidizing efficiency through 
direct effects on the ozone and OH photochemistry. Increases in water vapor concentration due to global 
warming will also lead to an increase in OH production. Observational data are totally inadequate for 
establishment of any changes in OH concentration and assessment of such changes must rely on model 
calculations. 

There have been a number of model evaluations of the impact of man's activities on tropospheric 
ozone and OH concentrations through increased emissions of methane, CO, NO x , and non-methane 
hydrocarbons. Initially, models with limited spatial resolution were employed, but increasingly one-, two-, 
and three-dimensional models have been developed to study both the chemical and dynamical aspects of 
tropospheric ozone (Crutzen, 1974; Fishman and Crutzen, 1978; Derwent and Curtis, 1977; Isaksen, 1981; 
Levy, Mahlman and Moxim, 1980; Logan et aL, 1981; Crutzen and Bruhl, 1989). Attention given to 
tropospheric modeling in previous ozone layer depletion evaluations has not been as detailed as that given 
to stratospheric modeling. Nevertheless, considerable progress has been made in theoretical studies towards 
resolving many of the issues relevant to an understanding of tropospheric sinks for alternative halocarbons 
and increased tropospheric ozone concentrations. 

The major conclusions of these studies can be summarized as follows: 

• The increased availability of NO x due to man's activities can change the Northern Hemisphere 
photochemistry from a net sink for ozone into a net source. 

• The Northern Hemisphere ozone budget is now dominated by man-made sources, particularly in 
the lower troposphere. 

• There may well have been a small decrease in the mean tropospheric OH concentration in the 
Southern Hemisphere since the pre-industrial era. The magnitude and direction of any change in 
the Northern Hemisphere is not clear. 

The accuracy with which any of these statements based on tropospheric modeling reflect what has 
actually happened in the real world depends on the adequacy and completeness of the model formulations 
together with their input assumptions. Significant progress in basic understanding has been achieved 
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towards the goal of predicting the extent and direction of man's influence on the oxidizing capacity of the 
troposphere. 


However, significant problems still remain. Potential sources of error or uncertainty in current one- 
and two-dimensional models arise principally through the following inadequacies of formulation: 

• the non-linear relationship between ozone production and precursor molecule concentrations, 
which creates difficulties in spatial averaging; 

• the short lifetime of NO* and its close coupling to a surface wet and dry deposition sink, which 
creates difficulties with large-scale transport modeling; 

• uncertainties in the latitudinal and seasonal distributions of the source strengths of methane, CO, 
NO x , and non-methane hydrocarbons; 

• uncertainties in the representation of tropospheric chemistry, including the representation of 
peroxy radical reactions, photodissociation processes (including the representation of clouds and 
aerosol scattering), nighttime chemistry and heterogeneous chemistry (including the removal of 
peroxy radicals on clouds, droplets and aerosols); and 

• inadequate observational data to test the models. 


In particular, the consideration of the first two areas listed above suggests that current models may over- 
estimate tropospheric ozone production from man-made sources. On the other hand, some modeling studies 
have demonstrated an important role for convective lifting of NO* from the boundary layer to the free 
troposphere with increased O^ production. Despite these considerable uncertainties, current models are 
able to account for the observed trends in tropospheric ozone concentrations in the Northern Hemisphere 
from pre-industrial times through to the present day. This requires a simultaneous and consistent treatment 
of man’s influence on the life cycles of methane, CO, NO*, and the non-methane hydrocarbons. Many of 
the relevant factors are poorly understood. However, there is no reason to anticipate that man’s influence 
on these trace gas life cycles will diminish over the next decade or so. It is pertinent to extrapolate current 
trends of these trace gas concentrations into the near future and to explore with the current models the 
influence of the changes on the oxidizing capacity of the troposphere in terms of OH and ozone. 

As worldwide industrial development continues in mid-latitudes and spreads to lower latitudes of the 
Northern Hemisphere, ozone concentrations are anticipated to grow throughout the Northern Hemisphere 
as a result of increased emissions of the precursors (Crutzen, 1988). The magnitude of any change in OH 
or ozone will be strongly scenario dependent and will be related directly to the pattern of change in the 
methane, CO, NO x , and non-methane hydrocarbon source terms (Isaksen and Hov, 1987; Thompson et 
al., 1989). Based on current perceptions of trends in trace gas concentrations, tropospheric ozone could 
increase by as much as 50% and the tropospheric mean OH concentration could decrease by as much as 
25% by the middle of the next century. 

The evaluation of future changes in tropospheric ozone and OH with current one- and two-dimensional 
models are subject to similar errors and uncertainties as those described above in the context of changes 
from the pre-industrial to the present-day atmosphere. Again, the non-linearities involved in relating OH 
and ozone to NO* concentrations give rise to concern. Models incorporating a more realistic description 
of the transport processes may indicate less sensitivity of tropospheric ozone and OH to changes in trace 
gas concentrations. 
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So far attention has been given to the direct influences on the oxidizing capacity of the troposphere 
and they can only be crudely represented in current one- and two-dimensional models. Several potentially 
important indirect influences have been proposed. These include: 

• increased UV-B penetration due to ozone depletion will increase tropospheric ozone photolysis 
rates, and consequently increase the mean OH source term (Liu and Trainer, 1987); 

• increased tropospheric temperatures may decompose methane clathrates and stimulate tundra 
methane emissions and biogenic hydrocarbon emissions (Ehhalt, 1988); and 

• increased water vapor mixing ratios may increase tropospheric OH concentrations. 

Understanding is growing steadily such that it should be possible in the near future to generate some 
consistent trace gas scenarios to aid in evaluating the influence of mankind on future tropospheric ozone 
and OH concentrations. 


4.3 OZONE DEPLETION POTENTIALS 

4.3.1 Background 

Recent consideration of international regulatory actions on the production of chlorofluorocarbons 
(CFCs) and other halogenated species has prompted significant interest in determination of the relative 
potential from such industrially produced compounds to affect stratospheric ozone and, more recently, 
global climate. The concept of relative Ozone Depletion Potentials (ODPs), introduced by Wuebbles ( 1981), 
has been adopted as a guideline or quick reference for estimating the relative potential for CFCs and other 
halocarbons to destroy stratospheric ozone. Several past papers have determined ODPs for selected 
chlorinated constituents (Wuebbles, 1983; Hammitt et al., 1987; Rognerud et al., 1989). This concept plays 
an important role in the implementation of the regulatory policies for fully halogenated CFCs adopted in 
the Montreal Protocol (UNEP, 1987). 

Several atmospheric modeling groups have recently reevaluated the ODPs for the relevant CFCs, the 
potential replacements (particularly several HCFCsand HFCs),and other industrially produced chlorinated 
and brominated compounds of interest being considered. The 1-D and 2-D global atmospheric models used 
in these evaluations incorporate the most recently available laboratory data on the chemical kinetics and 
radiative parameters of the species being examined. The discussion that follows primarily focuses on the 
recent studies by Fisher et al. (1989a, b), Connell and Wuebbles (1989), and Gillotay et al. (1989). Most of 
the calculations presented here are also contained, generally in greater detail, in the Fisher et al. (1989b) 
paper published in the AFEAS Report (Vol. II Appendix to this report). The Fisher et al. (1989a, b) papers 
are the result of a cooperative study among four modeling groups: Atmospheric and Environmental 
Research, Inc. (AER); ET. du Pont de Nemours and Company (DuPont); Lawrence Livermore National 
Laboratory (LLNL); and the University of Oslo. The Belgian Institute for Space Aeronomy (BISA) also 
contributed some 1-D model results to this report; these results were done separately and have not been 
as carefully compared as were the results of the other four groups. 

4.3.2 Definition of ODP 

Ozone Depletion Potential (ODP) has traditionally been defined as the steady-state ozone reduction 
calculated for each unit mass of a gas emitted per year (as a continuous release) into the atmosphere relative 


424 



OZONE DEPLETION POTENTIALS 


to that for a unit mass emission of CFC- 1 1 . This definition provides a single-valued relative measure of the 
maximum calculated effect on ozone of a given compound compared to the effect calculated for CFC-1 1 . 

Figure 4.3-1, taken from Fisher et al. (1989a), shows the resulting reduction in total ozone from the 
DuPont l-D model for 3 x HP kg of halocarbon emitted into the atmosphere during the first year, with no 
emission in subsequent years. As seen from the table insert in Figure 4.3-1, the time-integrated effects of 
the initial 5 x 1(F kg emissions agree closely to the normal steady-state determined ODPs. This confirms 
that ODP is a valid measure of the relative, cumulative effect on stratospheric ozone for each mass unit 
emitted into the atmosphere. 

It is also important to recognize the limitations in the ODP concept as defined above. As an example, 
ODP is defined in terms of the steady-state ozone change; it is, therefore, not representative of the relative 
transient effects expected for shorter-lived compounds during the early years of emission. During these 
early years, shorter-lived compounds may have changes in ozone much closer to the changes in ozone 
determined for longer-lived gases, such as CFC- 1 1, than would be expected based on their relative ODP 
values. 

Figure 4.3-2 shows the tropospheric concentration of two compounds following the onset of a constant 
emission level of each. One of the compounds has a 100-year lifetime (comparable to lifetimes of CFCs) 
and the other compound has a 5-year lifetime (comparable to lifetimes of hydrohalocarbons). We see that 
the ratio of the relative atmosphere concentration ratio does not approach steady-state until after about 



Figure 4.3-1. Calculated column ozone change in the DuPont 1-D model following a pulsed input 
(for one year) of 5.0 x 10 9 kg (Fisher et al., 1989 b) of specified gas. 
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Figure 4.3-2. Atmospheric concentrations of halocarbons with 5- and 100-year lifetimes following onset 
of a constant emission of each compound. 


400 years. Since the ODPs are based on steady-state relative effects which are (nearly) proportional to 
CFC concentrations, the calculated relative effects for short-lived compounds during the approach to 
steady-state can be larger than indicated by the ODP value. This is simply because the effect from the 
short-lived species has reached its full strength, whereas the effect from the long-lived species is still short 
of its steady-state value. 

Other uncertainties and limitations associated with the concept of and determination of Ozone 
Depletion Potentials will be discussed later. Included is further evaluation of time-dependent relative ozone 
depletion. 

In order to make the ODP definition consistent among models and at the same time have a conservative 
estimate of relative effects, the following criteria were selected for the model calculations included: 

1 . Depletion level— the calculations are based on emission rates of each compound required to give 
a modeled ozone depletion of approximately 1%. This value of depletion was selected in order 
to yield. results large enough to avoid the noise level inherent with numerical models, yet small 
enough to remain in the linear perturbation region. 

2. Trace gas levels — changing concentrations of other trace gases affect calculated tropospheric 
OH levels (Sze [ 1977] and Chameides et al. [1977]) and future depletions. However, we chose to 
base the calculations on current levels of C0 2 , CH 4 , CO, and N 2 0 due to the uncertainties in 
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future concentrations. The constant concentration assumption was chosen for clarity and sim- 
plicity. 

3. Chlorine levels — since the long-lived CFCs are present in today’s atmosphere and will affect 
chlorine chemistry over the time scales that hydrohalocarbons might be used, background 
halocarbon concentrations were assumed constant at current levels (yielding an asymptotic 
chlorine concentration of ~3 ppbv, with the exception of U. Oslo, who used 5.2 ppbv). 

4. Bromine chemistry — current levels of bromine compounds were included in the model chemistry 
where appropriate. 

Using the above provisions, the ODP for species a was then calculated as: 

ODP (a) = 80, (x) /80, (CPC-I I) 

where 80, (a) is the change in total ozone at steady-state per unit mass emission rate. CFC-I I is used as 
the reference gas. 

4.3.3 Model-Calculated ODPs 

Each compound considered enters the atmosphere at ground level and, in general, is removed by a 
combination of chemical processes: reaction with hydroxyl (OH) in the troposphere and the stratosphere, 
reaction with excited-state oxygen in the stratosphere, and photolytic breakdown by ultraviolet light in the 
stratosphere. Laboratory measurements indicate that CFCs (c.g., CFC-11, -12, -113, -114, and -115) are 
primarily destroyed by ultraviolet light in the stratosphere and, to a lesser degiee, by reaction with excited 
atomic oxygen; reaction with hydroxyl radical for these gases is unimportant. Hydrohalocarbons are 
removed by all three processes, but the predominant mechanism is reaction with OH in the troposphere. 
The reaction of the hydrohalocarbons with hydroxyl radicals leads to appreciably shorter atmospheric 
lifetimes compared to the CFCs. 


Reaction rate constants and photolysis cross sections required as model input were obtained consis- 
tently from a number of sources for most of the calculations presented here. Whenever possible, recom- 
mendations from the most recent evaluation by the NASA panel (DeMore et al., 1987) were used. Rate 
data for several key reactions have been reevaluated by Hampson et al. (1989), and Molina (1989), these 
were used when available. Some reaction rate parameters were obtained from the open literature (Davidson 
et al. 11978], Hubrich and Stuhl [1980]), while some data were only available from unpublished sources 
(Magid, private communication, 1988). Finally, tor some HFCs, no photolysis measurements could be 
found; since photolysis is of little importance for these species, neglecting this destruction mechanism has 
little impact on the derived lifetimes. Table 4.3-1 lists references and, when possible, the assumed reaction 
rates for destruction of each CFC, HCFC, and HFC. Several halons (namely, halon -1301, -1211, -1202, 
and -2402) are also considered in this report. The photolysis data needed for the halon species are based 
primarily on Molina et al. (1982). All models utilize their “best” representations for transport and radiation 
processes. 

Both one- and two-dimensional models have been used for these evaluations. One-dimensional models 
calculate the global-averaged altitudinal variation of the relevant atmospheric chemical processes. These 
models all share the advantage that they are easy to operate and evaluate. Two-dimensional (2-D) models 
allow examination of calculated effects over the full range of latitudes and seasons. As such, their added 
complexity is offset by better representation of nonlinear characteristics of global stratospheric transport 
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Table 4.3-1. Chemical rate data used in ODP calculations 


Species 

Photolysis 

Ref. 

OH Reactions 


O('D) Reactions 


Reaction Rate* 
A E/R 

Ref. 

Reaction 

Rate b 

Branching 

Ratio' 

Ref. 

CFC-11 

a 

— 

— 

_ 

2.3(- 10) 

0.75 

a 

CFC-12 

a 

— 

— 

— 

1.4(- 10) 

0.86 

a 

CFC-113 

b 

— 

— 

— 

2.0(- 10) 

0.80 

c 

CFC-II4 

b 

— 

— 

— 

1 6( - 10) 

0.80 

c 

CFC-115 

b 

— 

— 

— 

8.9( — 1 1) 

0.80 

c 

HCFC-22 

a 

I.2(- 12) 

1650 

f 

1 0( - 10) 

1.0 

f 

HCFC- 123 

d 

6.4( — 13) 

850 

f 

2.3( - 10) 

1.0 

f 

HCFC- 124 

d 

6.6(- 13) 

1250 

f 

l.0(- 10) 

1.0 

f 

HFC- 125 

— 

3.8( - 13) 

1500 

f 

0.5(-10) 

1.0 

f 

HFC- 134a 

e 

l.7<— 12) 

1750 

f 

0.5(-10) 

1.0 

f 

HCFC- 14 1 b 

d 

2.7(- 13) 

1050 

f 

1.51-10) 

1.0 

f 

HCFC- 142b 

d 

9.6( — 13) 

1650 

f 

1.41-10) 

1.0 

f 

HFC- 143a 

— 

2.6( — 13) 

1500 

f 

0.61-10) 

1.0 

f 

HFC- 152a 

e 

l.5( — 12) 

1100 

f 

1.01-10) 

1.0 

f 

CO, 

a 

— 

— 

— 

3.31-10) 

0.86 

a 

CH,CC1, 

a 

5.0( — 12) 

1800 

a 

3.181-10) 

0.80 

c 


^Reaction rate constant of form: k - A exp[ - E/(RT)], where A and k have the units of cm 3 /(molecule sec). 
h Reaction rate constant with units cm 3 /(molecu!e sec). 

‘'Fraction of 0( 1 D) disappearances proceeding through reaction channel, remainder pass through quenching channel to 0( 3 P) 
and no reaction with the halocarbon. 

Reference Key: 
a DeMore et al. (1987) 

b Hubrich and Stuhl (1980) 

c Davidson et al. (1978) 

d Molina (1989) 

e Magid (private communication, 1988) 

f Hampson, Kurylo, and Sander (1989) 


and radiative processes. Ozone changes calculated by 2-D models are averaged with respect to both latitude 
and seasons before calculation of ODPs. 


The calculated atmospheric lifetime of each gas species provides a key comparison of modeled results. 
Table 4.3-2 shows a comparison of lifetime values (of the CFCs, HCFCs, and HFCs) calculated by the 
models. In general, the calculated lifetimes for CFC-11 and CFC-12 are consistent with the lifetimes 
determined from the ALE-GAGE observations (Cunnold et al., 1986; Watson et al., 1988). However, the 
CFC-1 1 lifetimes from the 2-D models are on the low side of the uncertainty range, while 1-D results for 
the CFC-12 lifetime tend to be on the high side. The calculated HCFC and HFC lifetimes from all of the 
models are within the uncertainty range of those reported in section 4.2 as determined by Prather (1989) 
and Derwent and Volz-Thomas (1989) in the AFEAS Report (see Vol. II, Appendix). 


The Ozone Depletion Potentials calculated by eight models (from five research groups) for each of 
these species is given in Table 4.3-3. For each of the hydrogen-containing compounds, the calculated ODPs 
have been scaled by the ratio of the inferred lifetime for methyl chloroform (—6.3 years; see section 4.2) 
to the lifetime derived by each model. This scaling is an attempt to indirectly account for the different 
global tropospheric amounts of OH between models. The reaction with OH is a primary sink for the 
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Table 4.3-2. Atmospheric lifetimes (in years) calculated with one-dimensional and two-dimensional 


models 3 


Species 

Chemical 

Formula 


1-D Model Results 



2-D Model Results 

LLNL 

AER 

DuPont 

BISA 

Oslo 

LLNL 

AER 

DuPont 

CFC-ll 

CFCI, 

80.0 

60.0 

71.0 

85 

60.0 

52.0 

47.0 

46.0 

CFC-12 

CF 2 CI 2 

154.0 

125.0 

154.0 

166 

105.0 

101.0 

95.0 

118.0 

CFC-113 

CFCI 2 CF 2 CI 

96.0 

96.0 

117.0 

113 

101.0 

79.0 

— 

— 

CFC-114 

CF 2 CICF 2 CI 

209.0 

260.0 

319.0 


236.0 

197.0 

— 

— 

CFC-115 

CFjCICF, 

680.0 

690.0 

548.0 


522.0 

393.0 

399.0 

— 

HCFC-22 

CHCIF, 

20.0 

20.0 

16.0 


17.2 

15.0 

24.0 

12.7 

HCFC- 123 

CF,CHC1 2 

1.9 

2.1 

1.6 


1.7 

1.5 

2.4 

1.2 

HCFC- 124 

CF,CHCIF 

8.4 

8.8 

6.9 


7.4 

6.4 

10.0 

5.3 

HFC-125 

cf,chf 2 

37.0 

37.0 

25.0 


— 

27.0 

43.0 

19.0 

HFC-134a 

cf,ch 2 f 

21.0 

21.0 

16.0 


— 

15.0 

24.0 

12.5 

HCFC- 141b 

CCUFCH, 

8.9 

9.4 

7.8 


8.0 

6.9 

11.0 

5.8 

HCFC- 142b 

CCIFjCH, 

25.0 

25.0 

19.0 


20.9 

19.0 

28.0 

18.1 

HFC-143a 

CFjCHj 

54.0 

52.0 

42.0 


— 

40.0 

— 

— 

HFC-152a 

chf 2 ch, 

2.1 

2.3 

1.7 


— 

1.5 

2.7 

1.3 

Carbon Tetrachloride 

ecu 

73.0 

53.0 

61.0 

64 

52.2 

47.0 

40.0 

40.0 

Methvi Chloroform 

CH,CC1, 

7.4 

7.4 

6.0 


6.3 

5.8 

7.8 

4.7 

a The results of the LLNL, AER, DuPont, and Oslo groups have been published by Fisher et al. (1989a, b). 


Table 4.3-3. Ozone Depletion Potentials (ODP) calculated with one-dimensional and two-dimen- 
sional models, assuming scaling for HCFC ODPs by methyl chloroform inferred lifetime 3 


Species 

Chemical 

Formula 


1 -D Model Results 



2-D Model Results 

LLNL 

AER 

DuPont 

BISA 

Oslo 

LLNL 

AER 

DuPont 

CFC-ll 

CFC1, 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

CFC-12 

CF 2 C1 2 

1.0 

0.9 

1.0 

0.9 

0.9 

0.9 

0.9 

0.9 

CFC-113 

cfci 2 cf 2 ci 

0.8 

0.8 

0.9 

0.9 

0.9 

0.8 

— 

— 

CFC-114 

cf 2 cicf 2 ci 

0.8 

0.6 

0.8 

0.8 

0.6 

0.6 

— 

— 

CFC-115 

CFjCICF, 

0.4 

0.4 

0.5 

0.4 

0.4 

0.3 

0.4 

— 

HCFC-22 

CHC1F: 

0.05 

0.05 

0.04 

0.05 

0.04 

0.05 

0.06 

0.04 

HCFC- 123 

CF,CHC1 2 

0.13 

0.16 

0.13 

0.0l b 

0.013 

0.017 

0.022 

0.017 

HCFC-124 

CFjCHCIF 

0.016 

0.018 

0.017 

— 

0.018 

0.019 

0.024 

0.017 

HFC-125 

cf,chf 2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

HFC-I34a 

cf,ch 2 f 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

HCFC-14lb 

CCUFCH, 

0.07 

0.08 

0.07 

0.1 h 

0.09 

0.09 

0.11 

0.09 

HCFC- 142b 

CCIF.CH, 

0.05 

0.05 

0.05 

0.06 b 

0.06 

0.05 

0.06 

0.05 

HFC-143a 

CF,CH, 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

HFC-152a 

CHF 2 CH, 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Carbon Tetrachloride 

ecu 

1.1 

1.2 

1.1 

1.0 

1.2 

1.1 

1.0 

1.2 

Methvi Chloroform 

CHCCU 

0.10 

0.12 

0.10 

— 

0.14 

0.13 

0.16 

0.15 

Scaling assumes CH 3 CCI 3 lifetime of 6.3 years, based on discussion in Section 4.2. 

b Based on Gillotay et al. (1989). 
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hydrogen-containing species; differences in calculated tropospheric OH concentrations account for much 
of the model differences in chemical lifetimes for these gases. Therefore, amounts of OH account for much 
of the variation in ODPs derived between models. However, as noted by Fisher et al. (1989a, b), using 
globally averaged OH deduced from CH 3 CCL to calculate lifetimes for other HCFCs may not necessarily 
be accurate, particularly for those HCFCs with lifetimes substantially different from methyl chloroform. 

As seen in Table 4.3-3, the ODPs for each hydrogen-containing species, including the suggested 
replacement compounds, are much smaller than those for the CFCs. This difference largely reflects the 
reactivity of these species with OH, resulting in shorter chemical lifetime for the HCFCs compared to the 
CFCs. Several compounds, such as HFC-I34a and HCF-I52a, have ODPs of zero because they do not 
contain chlorine. 

In general, the ODPs derived using the different models compare extremely well with each other. 
However, significant differences between modeled ODPs still exist. The variation between models in the 
ODP for CH3CCI3, for example, is likely related to differences in the calculated distribution of OH with 
altitude and latitude between models. Although the two-dimensional models generally have a sounder 
physical basis for their formulation than existing one-dimensional models, there are no significant differ- 
ences between the l-D and 2-D model results. 

Table 4.3-4 gives the ODPs determined by three groups (AER and LLNL l-D models and U. Oslo 
2-D model) for several brominated halocarbons. Although the ODPs for these compounds are shown 
relative to CFC-ll for historical reasons, the ODPs for bromine-containing compounds should in reality 
be compared relative to each other because of the strong dependence to background chlorine levels in 
determining bromine effects on ozone. Therefore, in Table 4.3-4, Bromine Ozone Depletion Potentials 
(BODPs) are used for relative comparisons with halon-1301, the compound with the longest lifetime and 
largest ODP, as the reference. The major ozone destruction due to the halons comes from coupled chlorine- 
bromine chemistry in the lower stratosphere. Because of interactions between chlorine and bromine 
chemistry, the traditional ODPs are critically dependent on chlorine levels in the background atmosphere. 
Also, because current models underestimate the CIO concentrations in the lower stratosphere at high 
latitudes in early spring, the traditional ODP definition (relative to CFC-ll) may underestimate the true 
ozone depletion potential for these compounds. 


Table 4.3-4. Ozone depletions (ODPs) for brominated compounds as calculated in the AER and 
LLNL 1-D models and the U. Oslo 2-D model 



Atmospheric Lifetime 
(yrs) 


ODP" 



BODP B 


Common Name 

AER 

LLNL 

Oslo 

AER 

LLNL 

Oslo 

AER 

LLNL 

Oslo 

halon-1301 

81 

107.0 

72 

16.0 

13.2 

7.8 

1.0 

1.0 

1.0 

halon-121 1 

12 

15.0 

18 

3.5 

2.2 

3.0 

0.22 

0.17 

0.38 

halon-1202 


1.5 



0.3 



0.02 


halon-2402 


28.0 

23 


6.2 

5.0 


0.5 

0.64 


'"Relative to CFC-1 1, shown for historical purposes, Values will be underestimates if they account for effects at polar latitudes. 
Assumed upper stratospheric Cl* mixing ratio is 3 ppbv in the LLNL and AER models and 4 5 ppbv in the Oslo model. 
^Bromine Ozone Depletion Potentials (BODPs) defined relative to halon-1301, the longest brominated gas. 
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4.3.4 Uncertainties and Sensitivity Studies 


The above globally averaged ODP values are derived from model calculations based on present-day 
atmospheric conditions using current representations of chemical and transport processes. Although the 
ODPs agree reasonably well among the models, many uncertainties still exist. Establishing a strict criterion 
for estimating the overall uncertainty in the calculated ODPs is not a straightforward task. There are still 
many uncertainties associated with the treatment of atmospheric chemical, radiative, and dynamical 
processes in current models. Perhaps the single largest uncertainty is that none of the models used for 
calculating ODPs includes the chemical and dynamical processes causing the seasonal ozone losses asso- 
ciated with the ozone hole over Antarctica. The possible impact of this uncertainty is discussed below. 

Additional calculations have been carried out to test the sensitivity of the relative ODP effects to 
(I) levels of other trace gases, (2) seasonal and latitudinal dependencies and (3) assumed stratospheric 
transport processes. The relative time-dependent effects of different gases on total ozone is also examined. 


4.3.4.1 Ozone Hole Effects 

None of the ODP calculations consider the potential effect of heterogeneous chemistry in the lower 
stratosphere, particularly within the circulation vortex occurring at either pole during late winter and early 
springtime. Currently, inclusion of such effects in ODP calculations is premature since modeling of heter- 
ogeneous chemistry in general, and the polar phenomena in particular, are still in the early stages. 

Because of the apparent special chlorine processing and dynamics within the winter polar vortex, 
local Antarctic ODPs are expected to be larger than those shown in Table 4.3-3. Insofar as the observed 
long-lived tracer distributions, such as those for CFC-1 1 in the polar vortex, suggest that much of the total 
chlorine may be available there, then an upper limit on Antarctic ODPs can be determined by calculating 
the relative amounts of chlorine transported through the tropopause by the different gases. It can be argued, 
within existing uncertainties, that active chlorine at all altitudes within the stratosphere, and pet haps into 
the mesosphere, may be available for affecting ozone within the polar vortex. 

Values from four2-D models (AER, I.LNL, Oslo, DuPont) for the maximum relative chlorine loading 
into the stratosphere, termed Chlorine Loading Potentials (CLPs), are given in Table 4.3-5 based on the 
model calculated lifetimes. CLP. defined as the chlorine transported through the tropopause per mass 
emitted relative to CFC-1 1 , is proportional to the relative number of chlorine atoms per molecule and the 
atmospheric lifetime and is inversely proportional to the molecular weight. A set of reference lifetimes, 
based on the estimates used in scenario development in Chapter 3, section 3.2.1 tor the CECs and from 
the analysis in Section 4.2 for the HCFCs and HFCs, is given in Table 4.3-6 along with the corresponding 
CLP values. These values are as much as a factor of two to three larger than the derived ODP values. 

The reader may wonder why the values of the chlorine loading potential CLP are often larger than 
the ODPs for the same species. If all sources of stratospheric chlorine were equivalent in the photochemical 
reactivity sense once they reached the stratosphere, CLP values would be identical to the ODP values. 
Two major reasons account for the chlorine loading potential values being larger than the ODPs for many 
of the species. First, several of the CECs and HCFCs, such as HCFC-22 and HCFC-142b, are incompletely 
dissociated within the stratosphere. Second, differences in the dissociation reactivity produce different 
chlorine distribution for the various gases. 
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Table 4.3-5. Chlorine Loading Potentials (CLPs) from 2-D models for CFCs, HCFCs, and HFCs, 
scaled by lifetime of CH 3 CCI 3 ( = 6.3 yr) 


Species 

Oslo 

CLP 

LLNL 

CLP 

AER 

CLP 

DuPont 

CLP 

CFC-11 

1.0 

1.0 

1.0 

1.0 

CFC-12 

1.3 

1.5 

1.5 

1.9 

CFC-I13 

1.2 

1.1 



CFC-114 

2.1 

2.0 



CFC-115 

2.6 

2.2 

2.5 


HCFC-22 

0.15 

0.17 

0.22 

0.20 

HCFC-123 

0.017 

0.01 

0.03 

0.02 

HCFC-124 

0.04 

0.04 

0.06 

0.05 

HFC-125 

0.0 

0.0 

0.0 

0.0 

HFC-I34a 

0.0 

0.0 

0.0 

0.0 

HCFC-I4lb 

0.10 

0.11 

0.15 

0.13 

HCFC-I42b 

0.16 

0.18 

0.22 

0.20 

HFC-I43a 

0.0 

0.0 

0.0 

0.0 

HFC-I52A 

0.0 

0.0 

0.0 

0.0 

ecu 

1.0 

u 

1.0 

1.0 

CH.CCI, 

0.11 

0.12 

0.14 

0.14 


Table 4.3-6. Maximum relative Chlorine Loading Potential (CLP) for examined CFCs, HCFCs, HFCs, 
and other chlorinated halocarbons based on reference species lifetimes chosen to be 
compatible with available atmospheric measurements and modeling studies 


Reference" 


Species 

Lifetime (yrs) 

Chlorine Loading Potentials b 

CFC-II 

60.0 

1.0 

CFC-12 

120.0 

1.5 

CFC-II3 

90.0 

1.11 

CFC-114 

200.0 

1.8 

CFC-115 

400.0 

2.0 

HCFC-22 

15.3 

0.14 

HCFC-123 

1.6 

0.016 

HCFC-124 

6.6 

0.04 

HFC-125 

28.1 

0.0 

HFC-134a 

15.5 

0.0 

HCFC-14lb 

7.8 

0.10 

HFCF-142b 

19.1 

0.14 

HFC-I43a 

41.0 

0.0 

HFC-I52a 

1.7 

0.0 

CCL 

50.0 

1.0 

CH,CCI, 

6.3 

0.11 

^Lifetimes (e-folding time) are on estimates used in scenario development 
analysis in Section 4.2 for the HCFCs and HFCs. 

in Chapter 3, section 3.2.1 for the CFCs and from the 

b Chlorme Loading Potential is defined as the maximum chlorine transported across the tropopause per mass emitted relative to 

the same for CFC-1 1 ; it is proportional to 
to molecular weight. 

lifetime and the number of chlorine atoms per molecule, and is inversely proportional 
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At this time, until model studies of heterogeneous chemistry and validation of transport parameter- 
ization have been completed, it is difficult to put definite limits on the uncertainties of the model-calculated 
ODP values. The answer to this question depends on the extent of the influence that polar springtime ozone 
destruction will have on global ozone, and on the extent of the downward transport within the Antarctic 
and Arctic polar vortex. The preliminary results reported in Fisher et al. ( 1989b) indicate that comparison 
of model-simulated trace gas distributions with available observation may be useful to limit the uncertain- 
ties. The CLP values (i.e.. Table 4.3.6) provide an upper limit on ODP value within the Antarctic (and 
perhaps, the Arctic) region, while dilution effects, after breakup of the vortex, may extend higher ODP 
values to other latitudes. At this time, the ramifications of polar ozone depletion on the globally determined 
concept of ODPs are not clear. 

4.3.4.2 Sensitivity of ODPs to Other Trace Gas Levels 

ODP values have all been based on calculated ozone changes in a modeled atmosphere, assuming 
present-day levels of CFCs, CO;, CH 4 , CO, N,0, Br„, and other trace gases. Because it is important to 
quantify the effects of possible changes in the future trace gas concentrations on ODPs. we will present 
calculations of the sensitivity of the ODPs to trace gas concentrations that may occur within the next 
century if current concentration trends continue. 

Trace gas concentrations were varied from current concentrations as shown in Table 4.3-7. These 
concentrations were changed both individually and as an ensemble in the sensitivity calculations. Calculated 
changes in ozone for tested CFCs and HCFCs are relative to an atmosphere including the assumed trace 
gas perturbation. Two models were used to determine the sensitivity of atmospheric lifetimes and ODPs 
to trace gas concentrations. Calculations with the AER 1-D model assumed a background stratospheric 
chlorine content of 3 ppbv. Calculations with the Oslo 2-D model used a background chlorine amount of 

5.2 ppbv. Prior calculations with the Oslo model indicate little effect on the derived ODPs from the assumed 
chlorine background amount. 

Table 4.3-8 shows the calculated changes in atmospheric lifetimes of CFC-ll, CFC-12. HCFC-22, 
and HCFC-123 for the various sensitivity calculations. The largest effects are evident for the HCFCs in 
the cases where CH 4 and CO concentrations are perturbed. Both CH 4 and CO strongly influence tropo- 
spheric HO, chemistry; therefore the concentration of OH, and consequently the primary destruction of 
the HCFCs, is affected, resulting in the determined impact on atmospheric lifetimes. The relatively minor 
changes in lifetime of the CFCs arc primarily a result of increased photolysis rates for the trace gas 
perturbed atmospheres. 


Table 4.3-7. 

Change in ground level trace gas concentrations assumed 
in sensitivity studies 


Background 

Perturbed 

Trace Gas 

Atmosphere 

Atmosphere 

CO; 

340.0 ppmv 

680.0 ppmv 

ch 4 

1.6 ppmv 

3.2 ppmv 

CO 

100.0 ppbv 

200.0 ppbv 

N;0 

300.0 ppbv 

360.0 ppbv 

Br, 

20.0 ppt v 

40.0 pptv 
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Table 4.3-8. Sensitivity of calculated lifetimes to changes in trace gas values (based on AER 1-D 
results) 



CFC-1 1 

Lifetimes (years) 
CFC-12 HCFC-22 

HCFC-123 

Present-day atmosphere 

60 

125 

20 

2.1 

CO : (350 — ► 680 ppmv) 

61 

128 

20 

2.1 

CH, (1.6— » 3.2 ppmv) 

60 

126 

23 

2.5 

CO (100 ->200 ppbv) 

60 

125 

25 

2.7 

N,0 (300— » 360 ppbv) 

59 

123 

20 

2.1 

CH,Br (20 — » 40 pptv) 

59 

125 

20 

2.1 

All changes combined 

60 

127 

26 

2.9 


Table 4.3-9 presents ODPs resulting from these calculations. For both the AER 1-D model and the 
Oslo 2-D model, the largest sensitivity in the ODP for CFC-12 occurs for the assumed concentration change 
in N 2 0. As expected, ODPs for the HCFCs are affected most by change in CH„ and CO levels. 

The calculated sensitivities to background trace gas concentrations can be explained as follows. First, 
changes in atmospheric lifetimes will affect the ODP in that, as the lifetime increases/decreases, more/less 
of the CFC or HCFC is accumulated in the atmosphere at steady state for the same mass emitted. Secondly, 
the response of ozone to increases in chlorine can be different in the future atmosphere. However, this 
change in response will only affect the ODP estimates if the response of the individual CFCs or HCFCs is 
different from that of CFC-1 1 . An example of this involves the increased N 2 0 case: although the absolute 
changes in ozone for each of the CFCs and HCFCs were determined to be less than half of the ozone 
changes determined for the present atmosphere calculations, sensitivity of the derived ODPs is much 


Table 4.3-9. Sensitivity of ozone potentials to changes in trace gas values 



1 . 

AER 1-D Model 




CFC-12 

HCFC-22 

HCFC-123 


Present-day atmosphere 

0.92 

0.057 

0.019 


CO : (340— ► 680 ppmv) 

0.93 

0.051 

0.018 


CH 4 (1.6 — > 3.2 ppmv) 

0.93 

0.069 

0.023 


CO (100— > 200 ppbv) 

0.93 

0.070 

0.024 


N,0 (300 —* 360 ppbv) 

1.01 

0.071 

0.019 


CH 3 Br (20 — > 40 pptv) 

0.92 

0.055 

0.019 


All changes combined 

0.95 

0.076 

0.026 


2. 

Oslo 2-D Model (all calculations with 5.2 ppbv Cl y ) 



CFC-12 

HCFC-22 

HCFC-123 

HCFC-124 

Present-day atmosphere 

0.92 

0.046 

0.013 

0.018 

CO, (340 — » 600 ppmv) 

0.94 

0.049 

0.016 


CH, (1 .6 — * 3.2 ppmv) 

0.93 

0.058 

0.016 


N 2 0 (300 — » 360 ppbv) 

0.97 

0.042 

0.015 


Br, (18 — 33 pptv) 

0.91 

0.051 

0.016 

0.020 

All changes combined 

0.96 

0.062 

0.018 
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smaller. The increase in NO y will moderate the impact of increased chlorine by tying up a larger portion 
of the chlorine in the form of chlorine nitrate, C I N C) •, . This effect is most noticeable in the lower stratosphere 
where CINO, is more stable. The compensating effect should be more efficient for CFC-I I where more of 
the odd chlorine is found in the lower stratosphere. This accounts for the high ODP of CFC-22 in the 
increased N,0 case in the AER model; the Oslo 2-D model does not show this effect. 

In conclusion, both lifetimes and derived ODPs are moderately sensitive to significant changes in 
background trace gas concentrations. However, for the scenarios examined, the range in sensitivity of the 
ODPs is well within the overall uncertainty range for model-determined ODP values. Possible systematic 
errors, perhaps due to missing chemistry or other processes, may provide a much more significant impact 
on the derived ODP values. 

4.3.4.3 Variation in ODPs with Latitude and Season 

Another manifestation of the effect of differences between species on the chlorine distribution is the 
variations of the ODP values as a function of latitude and season. Several of the compounds show strong 
variations in their ODPs as a function ot latitude and season. 

Figure 4.3-3, based on calculations with the LLNL two-dimensional model (Connell and Wuebbles, 
1989; Fisher et al., 1989b), shows the calculated change in total ozone as a function of latitude and season 
for CFC-I 1. A constant emission rate of CFC-I I was assumed sufficient to give a steady-state global 



MONTH 

Figure 4.3-3. Calculated latitudinal and seasonal steady-state ozone change from emission of CFC-1 1 
necessary for global 1% change in total ozone (LLNL 2-D model). 
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decrease in total ozone of about one percent. As expected from previous 2-D analyses of CFC effects on 
ozone (e.g., WMO, 1985), the largest changes in ozone occur at the poles in late winter to early spring. 
Once the calculated changes are compared relative to CFC-11 and adjusted for differences in assumed 
emission rates, the resulting values indicate a generally variable dependency on latitude and generally 
weaker seasonal variations. In particular, the ODPs for CFC-12, HCFC-22, HCFC-124, and HCFC-142b 
have strong latitudinal variations, with the OOP increasing from lower tropical values to higher polar values 
by as much as a factor of three. The changes in total ozone at steady state from the LLNL model, and the 
derived ODPs from the LLNL, AER, Oslo, and DuPont 2-D models, for CFC-12 and HCFC-22 are shown 
in Figures 4.3-4 and 4.3-5, respectively. There are significant variations (more than a factor of two) with 
latitude and season in the ODPs derived for these two compounds. In both cases, the ODPs are smallest 
in the tropics and largest at the poles in late summer. The results from the models are in good general 
agreement regarding this pattern, although the LLNL model shows slightly more seasonal response for 
some species than the other models, particularly at high latitudes in the Southern Hemisphere. Other 
species, such as HCFC-123, HCFC-I4lb, and CH,CC1* show little variation with latitude or season (see 
Fisher et at. (1989b) and Connell and Wuebbles (1989) for figures and further discussion). 

The strong latitudinal ODP variation of some species (in contrast to the weak variation determined 
for others) warrants further discussion. The distinguishing difference between the species categories is 
based on the difference in the distribution of chlorine atoms which, in turn, depends on the altitudes for 
destruction of the respective species. Compounds that do not survive transport to the upper stratosphere 
for dissociation, and are therefore similar to CFC-1 1, show very little variation in their ODP with latitude 



MONTH 

Figure 4.3-4a. Calculated latitudinal and seasonal steady-state ozone change from emission of CFC-12 
necessary to give 1% change in total ozone (LLNL 2-D model). 
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MONTH 

Figure 4.3-4b. Calculated latitudinal and seasonal relative ozone depletion of CFC-1 2 (LLNL 2-D model). 



Figure 4.3-4c. Calculated latitudinal and seasonal relative ozone depletion of CFC-1 2 (AER 2-D model). 
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MONTH 

Figure 4.3-4d. Calculated latitudinal and seasonal relative ozone depletion of CFC-1 2 (Oslo 2-D model). 



Figure 4.3-4e. Calculated latitudinal and seasonal relative ozone depletion of CFC-1 2 (DuPont 2-D 
model). 
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Figure 4.3-5a. Calculated latitudinal and seasonal steady-state ozone change necessary to give 1% 
change in total ozone from emission of HCFC-22 (LLNL 2-D model). 



MONTH 

Figure 4.3-5b. Calculated latitudinal and seasonal relative ozone depletion of CFC-22 (LLNL 2-D model). 
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Figure 4.3-5c. Calculated latitudinal and seasonal relative ozone depletion from emission of HCFC-22 
(AER 2-D model). 



MONTH 

Figure 4.3-5d. Calculated latitudinal and seasonal relative ozone depletion from emission of HCFC-22 
(Oslo 2-D model). 
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Figure 4.3-5e. Calculated latitudinal and seasonal relative ozone depletion from emission of HCFC-22 
(DuPont 2-D model). 


and season. HCFC-123 is an example of such a species. On the other hand, species that persist into the 
upper stratosphere (even if the dominant removal is at lower altitudes) show strong latitudinal dependencies. 
CFC-12and HCFC-22 are two examples that show a significant latitudinal dependency. CFC-12is primarily 
dissociated at higher altitudes than CFC-11. HCFC-22, on the other hand, is primarily destroyed in the 
troposphere and lower stratosphere by reaction with OH radical, yet models and observations indicate that 
a significant fraction survives transport to the upper stratosphere. 

The resulting strong variations of ODPs with latitude occur in response to latitudinal differences in 
upper stratospheric chemistry on ozone destruction and the resultant effects of transport on both the CIO x 
produced and on the ozone being destroyed. The extent of these variations depends on the treatment and 
the strength of the modeled transport. It is not surprising, therefore, to find some differences in the ODP 
variations with latitude and season between models. The effects of heterogeneous chemistry and polar 
dynamical effects, not included in the calculations here, would lead to some modifications of these findings 
but should not qualitatively change the response. 

4.3.4.4 Sensitivity of ODP to Modeled Transport 

The question can next be raised as to the sensitivity of ODP values to model parameters — primarily 
transport processes. The treatment of transport processes is a source of significant uncertainty in 2-D 
models. Changing the model transport parameters would affect the ODP directly in three ways. First, the 
ODP values can change because of changes in lifetimes. The lifetime for CFC-1 1 is sensitive to the strength 
of the circulation in the stratosphere while the lifetime of HCFC-22, which is dominated by removal in the 
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troposphere, is less sensitive. Second, changing transport parameters would affect the distribution of 
chlorine atoms in the stratosphere. The responses of each species differ according to where the chlorine 
atoms are released. Finally, change in transport can also affect the response of ozone to chemical pertur- 
bations. 


Four sets of calculations were performed using different circulations and K yy (latitudinal eddy diffu- 
sivity) values to examine the sensitivity to variations in transport within the AER 2-D model. Calculated 
effects of CFC-1 1 and HCFC-22 emissions were determined in each case. These cases are designated by: 


Case 1 

standard 

standard (3 x 10 y cm 2 s ') 

Case 2 

standard 

small ( 1 x 10 K cm 2 s *) 

Case 3 

weak 

standard (3 x KTcnrs 1 ) 

Case 4 

weak 

small (1 x I0 X cm 2 s~ ') 


The weak circulation for cases 3 and 4 is chosen to test the response of ODP for HCFC-22 if a larger 
portion of the HCFC-22 is dissociated in the upper stratosphere. The circulation is adjusted so that the 
HCFC-22 molecules spend more time in the upper stratosphere, resulting in 90% dissociation of the 
molecules. Adopting smaller K yy values for cases 2 and 4 represents an attempt to obtain lower concentra- 
tions for species such as N.O, CFC-1 1, and CFC-1 2 in the polar lower stratosphere more in line with the 
observations obtained in the polar vortices during the Antarctic Airborne Ozone Experiment (AAOE) 
campaign. For simplicity, the smaller value was adopted globally and year round. Clearly, if the smaller 
values are restricted to the high latitudes and limited to certain seasons, the response could be somewhat 
different. 

Figure 4.3-6 shows calculated profiles of CFC- 1 1 from the four cases for August at 75°S compared to 
the AAOE measurements. The profiles at equinox for the Equator are shown in Figure 4.3-7. The cases 
with small K yy predict CFC-1 1 profiles that are in better agreement with observations at polar latitudes. At 
the same time, the concentration at the Equator calculated under case 2 is large compared to case 1. 
Corresponding graphs for HCFC-22 are shown in Figures 4.3-8. Figure 4.3-9 shows the calculated column 
abundances of O^ for the four cases. The results from cases 2 and 4 are distinctly different from observations 
with high columns predicted at high latitudes. Table 4.3-10 summarizes the calculated lifetimes and ODP 
values from the four cases. The calculated lifetime of CFC-1 1 is sensitive to the transport field, whereas 
the lifetime of HCFC-22 is less dependent. The ODP values for HCFC-22 for these cases reflect the impact 
of transport on the CFC-1 1 destruction process. Figure 4.3- 10a shows that while there is an increase in 
ODP values at high latitudes for case 2 relative to case 1 (Figure 4.3-5c), the decrease at the tropics more 
than compensates for it in the global ODP value. Corresponding results for cases 3 and 4 are shown in 
Figures 4.3- 10b and 4.3- 10c. 


In summary, the analysis presented here represents an attempt to quantify the sensitivity of the 
calculated ODP to model parameterization of transport. It is difficult to draw a definitive conclusion based 
on these limited number of simulations. The attempts to increase the ODP at high latitudes by using smaller 
K yy lead to much smaller ODPs at the tropics because of decreased eddy transport. It is possible that there 
may be other combinations of circulation and K yy that would change the ODP results significantly. These 
results did show that the same change in circulation would also have significant effects on the model- 
simulated ozone and other trace gas profiles. Comparison of these simulated results with available obser- 
vations would provide a way to validate the ODP values. 
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AUGUST 30 



Figure 4.3-6. Calculated CFC-1 1 profiles from transport sensitivity study at 75° south at equinox (AER 
2-D model). AAOE measurements indicated by cross hatch area. 


AUGUST 30 



Figure 4.3-7. Calculated CFC-1 1 profiles from transport sensitivity study at Equator at equinox 
(AER 2-D model). 


443 




OZONE DEPLETION POTENTIALS 


AUGUST 30 



MIXING RATIO 


Figure 4.3-8a. Calculated HCFC-22 profiles from transport sensitivity study at 75° south equinox (AER 
2-D model). 


AUGUST 30 



10 100 

MIXING RATIO 


Figure 4.3-8b. Calculated HCFC-22 profiles from transport sensitivity study at the Equator at equinox 
(AER 2-D model). 
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Figure 4.3-9a. i(a) Calculated total ozone abundance (Dobson Units) from transport sensitivity studies 
[case 1 of transport sensitivity study] (AER 2-D model), (b) Calculated total ozone abundance (Dobson 
Units) from transport sensitivity studies [case 2 of transport sensitivity study] (AER 2-D model). 




Figure 4.3.9b. ii(c) Calculated total ozone abundance (Dobson Units) from transport sensitivity stud- 
ies [case 3 of transport sensitivity studies] (AER 2-D model), (d) Calculated total ozone abundance 
(Dobson Units) from transport sensitivity studies [case 4 of transport sensitivity study] (AER 2-D 
model). 
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Table 4.3-10. Results of transport sensitivity study (AER 2-D Model) on ODP 
and CLP for HCFC-22 a 



Case 1 

Case 2 

Case 3 

Case 4 

CFC-ll lifetime 

47 years 

34 years 

51 years 

44 years 

HCFC-22 lifetime 

24 years 

22 years 

24 years 

23 years 

ODP 

0.07 

0.035 

0.055 

0.034 

CLP 

0.27 

0.35 

0.25 

0.28 

'These calculations do not include the CH 3 CCI 3 lifetime scaling included in Table 4.3-3. 



(a) MONTH 

Figure 4.3-1 0a. Calculated latitudinal and seasonal relative ozone depletion from emission of HCFC-22 [case 
2 of transport sensitivity study] (AER 2-D model). 
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(b) MONTH 

Figure 4.3-1 Ob. Calculated latitudinal and seasonal relative ozone depletion from emission of HCFC-22 [case 
3 of transport sensitivity study] (AER 2-D model). 



Figure 4.3-1 Oc. Calculated latitudinal and seasonal relative ozone depletion from emission of HCFC-22 [ case 
4 of the transport sensitivity study] (AER 2-D model). 
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4.3.4.S Time-Dependent Relative Ozone Depletion— Transient Effects 

One of the limitations associated with the ODP parameter is that it is based on steady-state pertur- 
bations; it does not consider the relative time-dependent effects of different constituents on ozone. While 
we have shown that ODP is equivalently a measure of the cumulative chronic effects of unit-emitted mass 
of a gas over its life cycle in the atmosphere, we are also interested in the transient response at a constant 
emission level. 

As mentioned earlier, ODP values determined at steady-state may not be representative of the relative 
transient effects expected for shorter-lived compounds during the early years of emission. During these 
early years, before one to two equivalent species lifetimes have passed, the ozone depletion per unit mass 
relative to CFC-11 can be much larger than the steady-state ODP value for some of the short-lived 
constituents. We note, however, that the changes in ozone occurring during the first few years after 
emissions are also relatively small compared to the steady-state ozone change. Relative time-dependent 
ozone depletions for several HCFCs and CFCs from the LLNL 1-D model are shown in Figure 4.3-1 1 and 
from the DuPont 1-D model in Figure 4.3-12. Results for the first few years are not shown because the 
derived changes in ozone were too small to be numerically accurate. 

The transient values for relative ozone depletion and the time it takes to approach the steady-state 
values depend on the time-dependent behavior of the calculated 0 3 decrease for the species examined and 
for the reference gas CFC-l 1. Since the atmospheric lifetimes of the HCFCs examined in this study are 
shorter than that of CFC-l 1 , the time constant for exponential decay to the ODP asymptote is determined 



Figure 4.3-11. Calculated time-dependent change in relative ozone column depletion following a step change 
in emission of the tested halocarbons (LLNL 1-D model). 
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Figure 4.3-12. Calculated time-dependent change in relative ozone column depletion following a step change 
in emission of the tested halocarbons (DuPont 1-D model). 


by the CFC- 1 1 lifetime. The magnitude of the transient relative ozone depletion depends on the stratospheric 
lifetimes and how rapidly the chlorine atoms are released after transport into the stratosphere. 

Several of the HCFCs, namely HCFC-123, HCFC-I4tb, and CH,CC1„ have much larger relative 
ozone depletions at early times compared to the ODP values. Each of these gases has relatively short 
reaction-time constants within the stratosphere — compared to CFC-11. As a result, chlorine atoms are 
quickly released by these gases and the full effect on ozone is reached much sooner than CFC- 1 1 . At longer 
times, the relative ozone depletions asymptotically approach the ODP value. 

Other species, such as HCFC-22, HCFC-124, and HCFC-142b, have long reactive time constants in 
the stratosphere. The time-dependent relative ozone depletions have much different behavior for these 
species. The relative ozone depletion builds gradually to a maximum value for these species, reflective of 
their long stratospheric time constants and the amount of time needed to build up stratospheric concentra- 
tions. After the first 20-40 years, the relative ozone depletion for these species gradually decays to the 
ODP value. 

For gases with longer lifetimes than the reference gas, no maximum in relative ozone depletion is 
noted. This can be seen by examining the results for CFC- 12 and CFC-1 15 in Figure 4.3-12. 

Figure 4.3-13 shows the relative time-dependent behavior of chlorine being transported through the 
tropopause for these same gases (taken relative to CFC-11). In all the cases, the time-dependent strato- 
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Figure 4.3-13. Calculated time-dependent relative chlorine loading following a step change in emission of the 
tested halocarbons (LLNL 1-D model). 


spheric chlorine loading is much larger initially compared with the eventual steady-state relative Chlorine 
Loading Potential (CLP) value (given previously in Table 4.3-6) if the lifetime of the species is less than 
the reference gas. 

In summary, time-dependent ozone depletion and chlorine loading values for the HCFCs are generally 
larger for the first 50 years than the ODP and CLP values defined at steady state. This, in part, arises from 
the longer atmospheric lifetime of the reference gas (CFC-1 1). Time-dependent relative ozone depletions 
for CFCs with lifetimes longer than CFC-I l show a monotonic increase to the steady-state ODP value. 

4.3.5 Summary 

One-dimensiona) and two-dimensional global atmospheric models have determined ODPs for a num- 
ber of halocarbons, including CFCs, other chlorinated compounds, several potential replacement hydro- 
halocarbons, and several brominated compounds. Although 2-D models have generally a sounder physical 
basis, there are no real differences between 1-D and 2-D results. In general, the ODPs for fully halogenated 
compounds, such as the CFCs, are much larger than those for the hydrogenated halocarbons, which include 
the potential replacement compounds considered. 

Brominated compounds should be compared with each other, because of the strong dependence of 
bromine effects on ozone to background chlorine levels. Bromine Ozone Depletion Potentials (BODPs) 
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are used for relative comparisons with Halon-1301, which has the longest lifetime and largest ODP, as the 
reference. 

Although the calculated ODPs agree reasonably well among models, many uncertainties still exist. 
None of the models used for calculating ODPs include the chemical and dynamical processes causing the 
seasonal ozone losses over Antarctica. Another uncertainty lies in the model-calculated OH, which is a 
major source of uncertainty for both lifetimes and ODPs of the HCFCs. 

Because of the apparent special chlorine processing and dynamics within the winter polar vortex, 
local Antarctic ODPs are expected to be larger than those derived. Insofar as the observed long-lived tracer 
distributions, such as CFC-I I in the polar vortex, suggest that much of the total chlorine may be available 
there, then an upper limit on Antarctic ODPs can be determined by calculating the relative amounts of 
chlorine transported through the tropopause by the different gases. These CLPs determined using assumed 
reference lifetimes (which generally agree with those in the models used here) can be as large as a factor 
of two to three times the derived ODP values. The ramifications of polar ozone depletion for global ODPs 
are not currently clear. 

The time-dependent relative ozone depletion values differ from the steady-state ODP values. The 
time-dependent values depend on the atmospheric lifetime and the transport time to the region of destruction 
of the gas. The shorter the stratospheric lifetime, the sooner the gas will impact stratospheric ozone and 
hence the higher the transient relative ozone depletion. For example, HCFC-123 has a shorter lifetime than 
CFC-11; its relative ozone depletion is largest soon after emission. Other gases in this category include 
HCFC-I41b and CH,CCI,. Species such as HCFC-22, HCFC-124, and HCFC-I42b have somewhat longer 
time constants in the stratosphere. Their relative ozone depletions build slowly to values (based on their 
time constants) as large as 0.2 and then decay slowly with time to the derived steady-state ODP value. 
Relative ozone depletion values for HCFCs are greater than ODP values even after 30 to 50 years. Time- 
dependent relative ozone depletions for CFCs with lifetimes longer than CFC- 1 1 show a monotonic increase 
to the ODP value. The ODP is the cumulative response; as discussed above, the transient response of 
relative ozone depletion may be larger than the ODP value at early times after emission. 

Several of the halocarbons indicate a strong latitude dependence in their ODP values and a generally 
weaker seasonal variation. In particular, ODPs for species such as CFC-12, HCFC-22, HCFC-124, and 
HCFC-I42b, which have stratospheric loss patterns very different from CFC-1 1 , produce strong latitudinal 
gradients in ODPs, with the largest ODPs near summer poles and smallest values in the tropics. The effects 
of heterogeneous chemistry and polar dynamical effects could modify these findings. 

Sensitivity analyses indicate that ODPs are affected to only a minor degree (< 20%) by assumed 
variations in background levels of NiO, CH 4 , COi, CO, total stratospheric chlorine, and total stratospheric 
bromine. 


4.4 HALOCARBON GLOBAL WARMING POTENTIALS 
4.4.1 Background 

The radiative perturbations and potential climatic effects due to the infrared absorption properties of 
chlorofluorocarbons (CFC-1 1 and CFC-12) have been assessed in WMO (1985). That assessment and those 
from accompanying studies (Wang and Molnar, 1985; Owens et al., 1985; Ramanathan et al., 1985, 1987; 
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Fisher et al., 1989c) demonstrated that the direct radiative forcing (i.e. , one in which all parameters are 
held fixed while the concentration of any CFCs is increased from 0 to 1 ppbv), consists of an increase in 
the downward longwave flux at the surface (surface radiative forcing) and a decrease in the net upward 
flux at the tropopause (surface-troposphere net radiative forcing). This imbalance results in a net increase 
in the amount of longwave radiation trapped within the Earth’s atmosphere system, thereby potentially 
enhancing the greenhouse effect. 


Since the 1985 reports, infrared spectral absorption data of several halogenated compounds have 
become available (Magid, private communication, 1988 and Gehring, private communication, 1988) which 
not only include CFCs, but also include HFCs and HCFCs. Interest in the latter two classes of compounds 
stems from their role as potential replacement compounds for CFCs. Laboratory measurements show that 
all the halogenated compounds possess significant absorption in the 500-1400 cm 1 (“window”) region. 
The absorption spectrum for any of these species typically is one of individual lines spaced very closely 
together in an arbitrary frequency interval giving the appearance of a continuum structure. The absorption 
cross sections of the HFCs and HCFCs also occur at similar frequencies in the “window” region and are 
comparable in magnitude to those for CFCs. 

Because the absorption bands of all replacement halogenated compounds fall in the infrared “window” 
region, the potential greenhouse effects need to be assessed. This need is emphasized by the fact that the 
clear portions of the Earth’s atmosphere have virtually no other constituents (excepting aerosols) with 
major absorption bands in the “window" region. In overcast conditions, clouds also have absorption bands 
in this region. Furthermore, the Earth’s outward-going radiation is at a maximum near the “window” 
region for the Earth s surface temperatures. Therefore, the presence of halogenated compounds could play 
an important part in determining the radiation balance in the clear portion of the atmosphere and a non- 
negligible role in the overcast portions. 


Besides the influence on the infrared radiative fluxes, there is yet another aspect of the perturbations 
by the halogenated compounds, as pointed out in WMO (1985). This is the heating perturbations induced 
in the upper troposphere and lower stratosphere which, owing to the large radiative time constant of those 
regions of the atmosphere, have the potential to increase the temperatures there. Any change in the 
temperatures near the tropopause region is of potential significance for the tropospheric-stratospheric 
exchange of water vapor. 


The effectiveness of greenhouse gases is accentuated by their lifetimes and in large part by their 
corresponding atmospheric abundances. 

In this report, we consider the most recently available laboratory spectral data of the halogenated 
compounds and assess the following: 

• the direct radiative forcing due to the halogenated compounds, as evaluated by both line-by-line 
(GFDL) and band models (AER, DuPont), 

• the changes in the surface temperature, as evaluated by the one-dimensional radiative-convective 
models (AER, DuPont), and 

• the relative changes in the surface temperatures taking into account the chemical lifetimes of the 
compounds expressed as Greenhouse Warming Potentials. 
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The GFDL line-by-line radiative transfer model has been developed from that described in Schwarz- 
kopf and Fels (1989). The band models employed in the AER and DuPont algorithms have been described 
earlier (Wang and Molnar, 1985; Owens et al., 1985). The absorption cross sections of the halogenated 
compounds are incorporated in a manner similar to the water vapor continuum. It should be noted that in 
the absence of a temperature-dependent dataset, none of these algorithms incorporate temperature depen- 
dence in calculated radiative effects. 

No attempt is made to calculate Greenhouse Warming Potentials for non-halocarbon gases such as 
carbon dioxide and methane. Due to their current atmospheric concentrations and the spectral locations 
of their infrared absorption bands, calculated global warming is not proportional to atmospheric concen- 
tration changes of these gases. In contrast, calculated warmings are linearly proportional to halocarbon 
concentrations. 

4.4.2 Definition of Halocarbon GWP 

Halocarbon Global Warming Potential (Halocarbon GWP) is defined parallel to the Ozone Depletion 
Potential. It is defined as the ratio of steady-state calculated warming for each unit mass of a gas emitted 
into the atmosphere (as indicated by the change in infrared flux at the tropopause) relative to the calculated 
warming for a mass unit of reference gas CFC-1 1. This definition was chosen as a representative measure 
of the potential of a compound to effect global warming, since: 

1 . It provides a measure of the cumulative effect on the radiative balance over its chemical lifetime 
for each unit released into the atmosphere (see below). 

2. The Halocarbon GWP yields a single value for each compound rather than a time-varying 
multitude of values. 

3. It provides a measure of the maximum calculated effect of a compound compared to the maximum 
calculated effect of CFC-1 1 on an equal mass basis. 

The first reason is perhaps the most important in that it estimates the cumulative effect on potential 
global warming of each unit released. An illustrative test was performed which quantified the chronic effect 
from a single pulsed release of test gas into the atmosphere. This particular test used the DuPont model to 
calculate cumulative global warmings over a 500-year time period following emissions of HCFC- 1 23, HCFC- 
22, and CFC-1 1. 

The calculated cumulative warmings are shown in Figure 4.4- 1 . For each case, the modeled warming 
is maximum very soon after the release and exponentially decays with a time constant equal to the 
atmospheric lifetime of the species. As seen in the insert table, the time-integrated warmings echo the 
relative values of the Halocarbon GWP calculated from steady-state figures. 

In order to make the definition of Halocarbon GWP consistent between models as well as a conser- 
vative estimate of relative effects, the following criteria have been selected: 

1. Trace gas levels — Changing the concentration of other trace gases could affect the calculated 
future equilibrium temperature rise from gases under evaluation here for two reasons. First, if 
there is overlap of absorption spectra, certain frequency bands would have less effect. Secondly, 
chemistry and therefore lifetime can be affected by perturbation of these chemicals. Current 
levels of C0 2 , CH 4 , N 2 0, 0 3 , and stratospheric H : Q were used in model calculations. 
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Figure 4.4-1. Calculated change in surface temperature following a pulsed emission of 5 x 10 9 kg of 
specified gas. 


2. Radiative forcing level — The calculations are based on gas concentrations that yield modeled 
responses large enough to avoid the “noise levels” of the numerical models, yet small enough 
to be in the linear response region. 

3. Reference gas — CFC-11 has been chosen as the reference compound for Halocarbon GWP 
calculations in order to have a reference material consistent for both Halocarbon GWP and ODP. 

Furthermore, additional parameters are useful in order to facilitate discussion: 

4. Net initial IR flux perturbation at the tropopause will be symbolically represented as dF (per unit 
mixing ratio in the troposphere). 

5. Specific Surface Temperature Change — The calculated surface temperature increase for a one 
part per billion (by volume) surface increase of any gas is defined as the specific surface temper- 
ature change, or symbolically dT„. 

6. The climate feedback factor (Dickinson, 1982), represented by X, is the ratio of the model- 
calculated surface temperature change to a perturbation in the net radiative forcing (discussion 
follows). 

The response of the Earth's surface atmosphere to a radiative forcing exerted by the presence of the 
halogenated compounds involves several feedback mechanisms and can be expressed as: 


dT s = X dF 
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where dT s is a measure of the climate response. The factor X depends on mechanisms not all of which are 
quantitatively understood. Nevertheless, as stated in WMO (1985). X can be expected to be in the range of 
0.25-1 K/W m 2 . Studies surveyed by Ramanathan et al. (1985) have two values ranging from 0.47-0.53, 
while the DuPont and AER models have a X of 0.64, and 0.81, respectively. It must be realized that the 
range in X stated above is not a rigid bracketing. Two important aspects of climate modeling severely 
restrict our present ability to better quantify X precisely. One is the lack of understanding of the response 
of the ocean circulation and the ensuing effects on the atmosphere. The second is the lack of knowledge 
of the changes in the cloud amounts and cloud properties due to changes in the halocarbon radiative forcing. 
Our knowledge in both the areas needs to improve substantially in order to better quantify X. 


The Halocarbon GWP definition resembles the ODP definition. For any gas, the general definition is: 


Halocarbon GWP = 


Calculated IR forcing due to compound X/ 
Emission rate (steady state) of compound X 
Calculated IR forcing due to CFC-1 1/ 
Emission rate (steady state) of CFC-1 1 


Note, since lifetimes are proportional to the ratio of atmospheric abundance to (molecular weight * 
emission rate), an equivalent form of this definition is: 


Halocarbon GWP = 


|dF(x) * Lifetime(x) / Molecular weight(x)] 
IdF(CFC-ll) * Lifetime(CFC-l 1) / Molecular weight(CFC-l 1)] 


(This formula has been used to calculate Halocarbon GWP from the line-by-line radiative model results.) 

Since radiative-convective models calculate a change in surface temperature, it is useful to express 
an equivalent form of the Halocarbon GWP, assuming the same value of X for all gases, as: 


Halocarbon GWP = 


|dTs(x) * Lifetime(x) / Molecular weight(x)] 

[dTs(CFC-l I) * Lifetime(CFC-l I) / Molecular weight(CFC-l I)] 


Note also, many of these gases have the potential of affecting heating rates indirectly as well, since they 
can chemically influence the distribution of ozone. Examination of model results indicates that this is a 
second order effect, at least two orders of magnitude below the IR effect and well below the sensitivity of 
these calculations (Wang et al., 1989). 


4.4.3 Model-Calculated Halocarbon GWPs 

A primary input to these radiative calculations are the altitudinal steady-state concentration profiles 
which are calculated using appropriate chemistry models. 

After the concentration profile is determined, the effect of each gas is calculated using a radiative 
model. This calculation accounts for the amount of energy absorbed by each IR gas over the radiation 
spectrum and accounts for the spectral overlap with other IR active gases. Quantitative infrared data for 
this input are available (see Fisher et al.. 1989c). 

Table 4.4-1 shows the net IR radiative flux at the tropopause ( = 12 km) for a 1 ppbv tropospheric 
mixing ratio, as calculated by each of the three models. Note that each of the model's calculations was 
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Table 4.4-1. 

Net IR radiative flux at the tropopause (at 12 km) [Wm 2 ] for 
1 ppbv tropospheric mixing ratio of each gas 

Species 

AER" 

DuPont b 

GFDL C 

CFC-1I 

0.175 

0.133 

0.35 (0.22) d , (0.30) e 

CFC-12 

0.248 

0.158 

0.42 

CFC-113 

0.223 

0.163 

0.42 

CFC-114 

0.258 

0.181 


CFC-115 

0.204 

0.164 

0.44 

HCFC-22 

0.151 

0.107 

0.29 

HCFC-I23 

0.140 

0.092 

0.28 

HCFC-124 

0.153 

0.108 


HFC-125 

0.189 

0.119 


HFC-I34a 

0.135 

0.095 

0.26 

HCFC-I4lb 

0.109 

0.076 

0.22 

HCFC-I42b 

0.144 

0.101 

0.27 

HFC-I43a 

0. 1 1 1 

0.087 


HFC-I52a 

0.092 

0.059 

0.17 

CCI 4 

0.080 

0.063 


CH,CCI, 

0.038 

0.033 


|2 x CCh] 

4.41 

3.87 

5.68 


Note difference in base assumptions: 

'AER calculations based on 
—global mean conditions 
— 48 5% cloud cover 

— cloud with spectral optical thickness at 0.55 m of 8.9. (cloud albedo = 0.54} 

‘ DuPont calculations based on 
—global mean conditions 
— 50% cloud cover 
— cloud albedo = 0.5 

■ GFDL calculations based on 
— mid-latitude, summer atmosphere 
— clear sky conditions 
— uniform gas mixing ratio 

’CFC-I 1 only in troposphere 

^CFC-1 1 only in troposphere, with 50% black cloud cover at —4 km 


based on different cloud assumptions (GFDL used a clear sky, DuPont used 50% cloud cover (albedo = 
0.5], and the AER was based on control climate cloud field with 48.5% cloud cover). Note also that the 
GFDL model assumed mid-latitude summer radiative conditions, whereas both AER and DuPont models 
are based on globally and annually averaged radiation. 

Specific surface warmings are calculated using the AER and DuPont radiative convective models and 
are shown in Table 4.4-2 for each of these gases. Again differences in the \ factors between these models 
yield a systematic difference in the response of surface warming to net radiative flux changes for each 
compound. Note that HCFCs and HFCs on the average have lower specific surface warming values than 
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Table 4.4-2. Specific surface temperature (dTs) increases and lambda (A) values resulting from 
1 ppbv of each gas 


Species 

AER 

dTs 

(°K/ppbv) 

Du Pont 

X 

(°K/ppbv/Wm 

AER 

-*> 

Du Pont 

CFC-11 

0.135 


0.084 

0.771 

0.632 

CFC-12 

0.202 


0.102 

0.815 

0.647 

CFC-113 

0.174 


0.103 

0.780 

0.632 

CFC-114 

0.208 


0.115 

0.806 

0.635 

CFC-11 5 

0.170 


0.107 

0.833 

0.652 

HCFC-22 

0.124 


0.070 

0.821 

0.650 

HCFC-123 

0.111 


0.059 

0.793 

0.644 

HCFC-124 

0.126 


0.070 

0.824 

0.645 

HFC-125 

0.160 


0.078 

0.847 

0.654 

HFC-I34a 

0.1 14 


0.061 

0.844 

0.647 

HCFC-I41b 

0.086 


0.048 

0.789 

0.637 

HCFC-142b 

0.120 


0.066 

0.833 

0.651 

HFC-143a 

0.092 


0.054 

0.829 

0.625 

HFC-152a 

0.076 


0.038 

0.826 

0.649 

CCI 4 

0.062 


0.040 

0.775 

0.628 

CH,CC1, 

0.027 


0.020 

0.710 

0.618 

[2 x CO,] 

[3.08J 


11.651] 




CFCs by approximately 40%. (Net dF values are also lower by -40%. ) Although total absorptances are 
about equivalent for each class of compounds, overlap with absorptance bands of other infrared active 
gases diminishes the efficiencies of these gases as greenhouse gases. Most of the hydrogenated halocarbons 
have bands that overlap the bands of both CH 4 and N 2 0 as well as water vapor above 1250 cm” 1 . 

Table 4.4-3 shows the Halocarbon GWP values calculated by each of the three models. Reference 
lifetimes are based on the estimates used in scenarios developed in Section 3.3.1 for CFCs and from the 
analysis in Section 4.2 for HCFCs and HFCs. 

Halocarbon GWPs for fully halogenated compounds are larger than those for the hydrogenated 
halocarbons. Fully halogenated CFCs have Halocarbon GWP values ranging from 1.0 to 7.5, whereas 
HCFCs and HFCs range from 0.02 to 0.7. 

Halocarbon GWP values differ between species due to differences in infrared absorbances and 
differences in lifetimes. The range of absorbances is approximately 4, while the lifetimes range by a factor 
of 250. Thus, the range of 400 among the Halocarbon GWP values is primarily a result of differences in 
lifetimes. Additionally, the Halocarbon GWP for halon-1301 has been calculated based on specific surface 
warmings reported in WMO 1985 and a relative lifetime (to CFC-11) of 1.3 (from LLNL l-D and Oslo 2-D 
models). This Halocarbon GWP value is 1.4. 
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Table 4.4-3. Halocarbon global warming potentials 


Species 

Reference* 
Lifetime (Yrs) 

AER b 

DuPont b 

GFDL b 

CFC- 11 

60.0 

1.0 

1.0 

1.0 

CFC- 1 2 

120.0 

3.4 

2.8 

2.7 

CFC-113 

90.0 

1.4 

1.4 

1.3 

CFC- 114 

200.0 

4.1 

3.7 


CFC- 115 

400.0 

7.5 

7.6 

7.5 

HCFC-22 

15.3 

0.37 

0.34 

0.34 

HCFC-123 

1 .6 

0.02 

0.017 

0.019 

HCFC-124 

6.6 

0.10 

0.092 


HFC-125 

28.1 

0.65 

0.51 


HFC-134a 

15.5 

0.29 

0.25 

0.26 

HCFC-I4!b 

7.8 

0.097 

0.087 

0.096 

HCFC-I42b 

19.1 

0.39 

0.34 

0.34 

HFC-143a 

41.0 

0.76 

0.72 


HFC-152a 

1.7 

0.033 

0.026 

0.029 

CC1 4 

50.0 

0.34 

0.35 


CH,CCIj 

6.3 

0.022 

0.026 


‘"Lifetimes are based on estimates used in scenario development (Chapter 3, Section 3.2, 1 ) for CFCs and from the analysis in 
Section 4.2 for HCFCs and HFCs. 

b AER and DuPont results are based on surface temperature perturbations calculated using radiative-convective models; GFDL 
results are based on line-by-line determined radiative forcing. 


4.4.3.1 Sensitivity to Trace Gas Levels 

The Halocarbon GWP values have all been based on calculations that assume present-day values for 
concentrations of trace gases of CCb, CH 4 , CO, N 2 0, and CFCs. Because the Halocarbon GWP values are 
for consideration with future atmospheres different from today, we have examined the sensitivity of these 
parameters to changes in C0 2 and CH 4 at levels that might be achieved if current trends continue. 

Trace gas changes and the resulting calculated changes in lifetimes, surface temperature rises, and 
Halocarbon GWP values (from the AER model) are shown in Table 4.4-4. Calculated changes in surface 
temperature for tested CFCs and HCFCs were compared to a reference atmosphere including the assumed 
gas perturbation (note the Halocarbon GWP values are different than listed above, since they are based on 
model-calculated lifetimes rather than reference lifetime values). 

As seen in Table 4.4-4, trace gas changes have little effect on Halocarbon GWP values. C0 2 has the 
greatest effect on CFC radiative forcing because the albedo feedback is weaker in the warmer atmosphere. 
As a result, the radiative forcing for CFC-1 1 and -12 is somewhat weakened. Methane affects the chemical 
lifetimes of the HCFCs, resulting in slightly greater Halocarbon GWP values for these compounds. Never- 
theless, these relatively major perturbations in trace gas levels have only minor effect on calculated 
Halocarbon GWP values, indicating that the GWP formalism is robust over the range of trace gases projected 
for the next century. 
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Table 4.4-4. Sensitivity study— variation in trace gas levels and impact on halocarbon GWP values 
(AER 1-D model) 



Atmosphere 

CFC-1 1 

Species 

CFC-1 2 HCFC-22 

HCFC-123 

Lifetime 

Present Day 

60 

125 

20 

2.1 

(yrs) 

CH 4 (1.6— » 3.2 ppbv) 

60 

126 

23 

2.5 


C0 2 (340 —* 680 ppm) 

61 

128 

20 

2.1 

dT s 

Present Day 

0.14 

0.20 

0.12 

0.11 

(°K/ppbv) 

CH 4 (1.6-* 3.2 ppbv) 

0.14 

0.20 

0.12 

0.11 


CO, (340 -» 680 ppm) 

0.12 

0.17 

0.12 

0.11 

GWP 

Present Day 

1.0 

3.5 

0.49 

0.026 


CH 4 (1.6— * 3.2 ppbv) 

1.0 

3.5 

0.56 

0.031 


C0 2 (340 -» 680 ppm) 

1.0 

3.4 

0.49 

0.027 


4.4.3.2 Time Dependency of Relative Global Warming 

Since the Halocarbon GWP parameter is based on steady-state effects, it does not describe the relative 
time-dependent effects of constituents on warming. Even though the Halocarbon GWP is an equivalent 
measure of the cumulative warming during its lifetime for each unit mass emitted, the transient response 
to a constant emission level is also of interest. 

The calculated warmings for a number of halocarbons are shown in Figure 4.4-2. As seen, the 
calculated warming reaches an asymptote rapidly for the HCFCs, but requires longer periods to approach 
steady-state for CFCs. These response patterns echo the respective patterns for increases in atmospheric 
abundances for each species, seen in Figure 4.4-3. 

Relative warmings are shown in Figure 4.4-4. For HCFCs, the relative effects are at a maximum at 
very short times. One can easily show (using L'Hospital’s rule) that the initial relative value is the ratio of 
the value of the specific surface temperature increase ([°K/ppbv] / molecular weight— relative to the same 
ratio for CFC-1 1). Thus initially, while the individual effects are small, the relative effects are on the order 
unity. However, as atmospheric concentrations build and chemistry differences affect the relative amount 
in the atmosphere, the relative effects either grow or decrease depending on whether the lifetimes are 
longer or shorter than that of CFC-1 1. As seen in Figure 4.4-4, the HCFCs have lifetimes shorter than the 
lifetime of the reference gas and have relative effects that asymptotically approach the Halocarbon GWP 
value with a time constant equal to the lifetime of CFC-1 1 . Longer-lived species have relative effects that 
grow with time, asymptotically approaching their Halocarbon GWP value with time constant for their own 
lifetime. 

4.4.4 Uncertainties 

Uncertainties in the effectiveness of gases to produce global warming fall into two general classes — 
those that are generalized for the entire greenhouse effect and those that are specific to individual species. 
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Figure 4.4-2. Change of calculated warming following a step change of emission of specified gas at 
5 x 10® kg/yr (DuPont 1-D model). 



Figure 4.4-3. Column of (hydro) chlorofluorocarbons following a step change of emission of specified 
gas at 5 x 10® kg/yr (DuPont 1-D model). 
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Figure 4.4-4. Calculated relative warmings following a step change of emission of specified gas [CFC-1 1 
reference] (DuPont 1-D model). 


There are a number of uncertainties in the modeling of greenhouse warming. These concern the 
radiative properties of the Earth’s surface-atmosphere system, such as changes in the surface and ice cover 
albedos, and changes in cloud cover and composition. Changes in the temperature structure of the atmo- 
sphere will affect the convective patterns and chemistry of the stratosphere. The coupling of oceans (as 
heat reservoirs) and ocean currents to surface temperature changes will also affect the timing and location 
of the warming. Research is being carried out to understand these questions, which apply to all gases that 
affect the future radiative balance of the Earth. 

Uncertainties also exist regarding the individual species and their influence on the radiative balance. 
Temperature dependence of the absorptances of individual species need to be resolved and modeled. The 
chemical processes affecting both lifetimes and atmospheric profiles are also the subject of continuing 
research. Resolution of questions related to the general greenhouse warming modeling will directly affect 
timing and magnitude of global warming, whereas resolution of the radiative and chemical parameterizations 
for halocarbons will have a direct impact on the Halocarbon GWP values for these species. 

4.4.5 Summary 

Halocarbon Global Warming Potentials have been defined and calculated in order to allow estimates 
of the relative environmental effects of halocarbons. The results presented here indicate that the Halocarbon 
GWPs depend primarily on the atmospheric lifetime of the compounds and to a lesser degree on the 
molecular IR absorption characteristics. 
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The reduction in Halocarbon GWP that might be expected due to use replacement of a CFC by a 
hydrohalocarbon can be estimated by taking the ratio of the Halocarbon GWP of the hydrohalocarbon to 
the Halocarbon GWP of the CFC it would replace. For example, the reduction in Halocarbon GWP in 
replacing uses of CFC-12 by HCFC-134a is (0.26 ± 0.030)/(3.07) = 0.085 ± 0.01. Of course, the relative 
quantities of the compound required in the use application must also be taken into account. 

Although the Halocarbon GWP values reported here agree reasonably well among models, uncer- 
tainties in the values still exist because of the uncertainties in modeled chemistry and dynamics and their 
direct effect on the chemical lifetimes of these compounds. 

Halocarbon GWP values appear to be reasonably robust parameters since their calculated values are 
nearly insensitive to assumed concentrations of other radiative gases. The minor shifting of the Halocarbon 
GWP values is primarily influenced by the changes in calculated lifetimes and therefore the abundance in 
the atmosphere. 

Calculated time-dependent relative global warmings for Halocarbons are initially on order unity but 
decrease or increase depending on whether their lifetimes are shorter or longer than that of the reference 
gas. At longer times, the relative global warmings asymptotically approach the Halocarbon GWP values. 
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